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How  to  U«e  thi»  Guide; 


1.  Get  e  printout  of  the  FCHITRAM  SOUECE  or  coapiled  vereione  of  the 
progreae  end  eubroutines  that  aeke  up  this  eyetea. 

2.  Ihie  User's  Guide  is  s  SUPPLEMENT  to  the  bulk  of  the  systea 
docuaentetion  included  as  extensive  coaaenting  of  the  progrsa  coding. 

3.  If  you  want  to  know  about  the  GENERAL  LAYOUT  of  the  systea,  you  can 
find  that  inforaation  conveniently  either  in  this  guide  or  by  looking  at 
the  aainlines  of  the  four  prograas  that  make  up  the  systea. 

4.  If  you  want  to  know  about  the  reconaended  LIMITATIONS  of  this  systea, 
this  docuaent  will  answer  those  kinds  of  questions. 

5.  If  you  want  to  know  WHY  A  COMPONENT  MAS  MODELLED  THE  HAY  IT  WAS,  this 
docuaent  should  answer  those  questions,  too. 

6.  This  guide  also  includes  some  discussion  of  the  GENERAL  CONVENTIONS  we 
followed  in  developing  this  systea. 

7.  You'll  have  to  look  at  the  PROGRAM  PRINTOUTS  to  find  the  DETAILED 
APPROACH  we  used  in  modelling  the  engine  coaponents  and  performance. 


Iv 


Chapter  I 

INTROOUCTIOIf:  What 'a  It  All  About? 


What  i»  thi«  «y«te«  suppo««d  to  do? 

Thia  Turbofan  Engine  Technology  Evaluation  Syatea  is  intended  to  be  a 
eoaputer  age  "back-of-the-envelope”  calculation  tool  for  use  in  evaluating 
the  relative  payoffs  of  coapeting  gas  turbine  engine  technologies,  d 
systea  capable  of  perforaing  all  the  tasks  required  to  fit  that 
deacription  aust  be  versatile  and  that's  one  thing  this  systea  is.  To  be 
useful,  we  thought  the  systea  should  also  be  easy  to  use  and  understand. 
The  page  titled  "How  to  Dae  this  Guide"  is  intended  to  give  you  a  head 
start  toward  using  and  understanding  everything  you  want  to  know*  about 
this  systea. 

Who  wrote  it? 


This  systea  was  developed  to  fulfill  one  of  the  goals  of  a  basic  research 
project  at  the  Frank  J.  Seiler  lesearch  Laboratory  (Air  Force  SystMs 
Coaaand),  a  tenant  organisation  located  at  the  United  States  Air  Force 
Aeadesqp  in  Colorado  Springs,  Colorado.  The  authors  included  two  full-tiae 
researchers,  Haj  Arthur  E.  Fanning,  and  Capt  Victor  J.  Villhard;  two 
reserve  officers,  Capt  Williaa  Oieterich  and  Capt  C.  Bruce  Haraon;  two 
college  students  who  worked  at  the  Isb  between  seswsters,  Mr.  Dwayne  Davis 
and  Mr.  A.  Scott  Mediavilla;  and  an  Air  Force  Academy  graduate  awaiting 
Pilot  Training,  2Lt  John  H.  Wood.  It  took  us  from  Fall  1981  to  Winter 
1983  to  completely  write  and  documnt  this  system. 

What  proarams  make  up  the  system? 

SETUP:  an  interactive  graphics  program  used  to  select  and  input 

design  cycle  parameters  for  the  engine  and  the  flight 
conditions  at  which  performance  is  to  be  calculated. 

TBBOFH:  the  turbofan  engine  performance  simulation  program 

ENCOM:  an  interactive  engine  vising  and  comparison  program  that  also 

allows  selection  of  performance  parameters  for  graphs 

GBAFH:  an  interactive  graphics  program  for  the  creation  of 

publication  quality  graphs  of  almmst  any  calculated 
parameter  versus  any  other. 

Are  these  programs  easy  to  understand  and  useT 

All  four  prograsw  in  this  system  are  modular.  That  is,  they  are  made  up 
of  calls  to  many,  many  short,  easy  to  uadarstaad  subroutines.  Bach 
sobroutina  is  co^letaly  internally  documented  and  eaplalned  in  the 
programming  coda.  Moat  are  also  mentioned  in  this  User  Guide. 

Bow  might  I  use  this  system? 

X  hope  you'll  find  thia  syetam  useful  in  your  tasks  related  to  turbofan 
deeign  and  perfnrnante  analyaie  as  well  as  in  training  now  people 
who  mi^t  bo  unfamiliar  vith  engine  performance  calculations  or  trenda. 
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What'«  in  thia  document  that  isn't  in  the  prograaminE  code? 

Thia  Uaer  Guide  deacribes: 

o  the  general  layout  of  the  syaten;  in  other  words,  the  picture  that 
covers  just  a  little  more  than  the  level  that  describes  the  purposes  and 
workings  of  the  main  programs, 

o  limitations  of  the  system  models  and  suggested  constraints, 

o  the  general  approach  used  in  calculating  gas  properties  throughout 
the  engine  simulation  model,  TRBOFM, 

o  why  the  engine  components  were  modeled  the  way  they  were  in  TRBOFN, 

o  what  output  products  you  can  get  from  this  system, 

o  how  the  performance  calculated  by  TRBOFN,  the  engine  simulation 
program  in  this  system,  compares  to  SMOTE,  an  accepted  standard  turbofan 
engine  simulation  program. 


GENERAL  INFORMATION  ON  THE  PROGRAMS; 

SETUP  is  an  interactive  graphics  program  which  allows  the  user  to  select 
and  change  engine  design  prameters  and  off-design  flight  conditions  for 
TRBOFN  to  use  in  calculating  engine  performance.  On  each  run  of  SETUP, 
you  can  design  up  to  five  engines  to  be  run  by  TRBOFN.  SETUP  creates  a 
data  file  that  TRBOFN  reads. 

TRBOFN  is  a  parametric  engine  simulation  program  capable  of  calculating 
performance  trends  for  both  mixed  and  separate  flow  engine  designs.  Gas 
properties  are  calculated  by  subroutines  which  reproduce  the  calculations 
Keenan,  Chao  and  Kaye  outline  in  their  book.  Gas  Tables.  TRBOFN  has  a 
control  schedule  built  into  it  to  simulate  performance  of  the  engine  at 
throttled  power  settings.  Calculated  perforaunce  consiata  of 
temperatures,  pressures,  Mach  numbers  and  gas  properties  at  about  50 
stations  through  the  engine,  engine  thermodynamic  efficiencies,  thrust  and 
specific  fuel  consumption. 

ENCOM  performs  sizing  and  comparison  functions  on  up  to  five  engines  whose 
performance  was  calculated  by  TRBOFN.  It  scales  the  physical  size  of  each 
engine  so  that  each  can  propel  the  same  airplane  through  a  set  of  up  to 
five  user-chosen  maneuver  conditions.  For  instance,  you  might  want  the 
airplane  to  be  able  to  (1)  sustain  a  6g  turn  at  Mach  0.9  at  20,000  ft 
altitude,  (2)  be  able  to  produce  a  Ig  acceleration  at  sea  level  for 
takeoff,  (3)  be  able  to  dash  at  Mach  2.3  at  10,000  ft  and  (4)  cruise  at 
Mach  0.85  at  40,000  ft.  ENCOM  includes  a  drag  polar  for  a  fighter  type 
airplane  which  can  easily  be  modified  to  fit  any  airplane  you  want.  It 
determines  the  thrust  requirements  for  each  of  the  maneuvers  you  specified 
from  the  drag  parameters.  It  then  scales  the  physical  size  of  each  of  the 
engines  to  produce  enough  thrust  at  each  flight  condition  to  meet  each 
maneuver  requirement.  Vftiichever  maneuver  required  the  largest  size  for 
each  of  the  five  engines  is  the  ’’critical  maneuver  requirement"  for  that 
engine.  The  relative  size  information  for  each  of  the  five  engines  is 
then  stored  and  output.  ENCOM  can  be  used  to  determine  the  optimum  engine 
cycle  in  terms  of  physical  size  for  a  given  mission  by  specifying  certain 
oianeuver  conditions  and  looking  at  the  relative  sizes  of  engines  with  a 
range  of  cycle  parameters.  But  it  even  goes  a  step  further  by  allowing 
the  user  to  compare  the  relative  performance  of  the  engines  flying  through 
an  eleven-leg  mission.  Using  this  section  of  the  program,  you  can 
optimize  the  engine  based  on  TSFC  or  thrust  margins  available  or  any  other 
criterion  you'd  like  to  use.  ENCOM  also  collects  parameters  from  the 
TRBOFN-generated  engine  performance  data  files  for  plotting  so  you  can 
examine  trends  and  determine  optimum  cycle  designs  using  graphic  displays 
to  present  complex  trends  and  relations  in  an  easily  understandable  and 
readily  usable  form. 

GRAPH  is  the  program  that  allows  you  to  plot  performance  trends,  etc.  It 
is  an  interactive  graphing  program  which  uses  Tektronix'  Interactive 
Graphics  Library  (IGL).  The  graphing  program  also  depends  on  features 
unique  to  Tektronix  411X  series  terminals.  GRAPH  could  easily  be  modified 
to  work  on  almost  any  other  Tektronix  terminal.  It  would  probably  be  less 
versatile  but  IGL  does  include  subroutines  to  allow  other  terminals  to 
esailate  auny  ot  the  411X  features.  Obviously,  if  you  don't  have  access  to 
ICL  and  you  do  have  access  to  another  graphics  package  you  could  write  a 
program  using  it  to  replace  GRAPH.  The  only  constraint  is  that  it  would 
have  to  be  written  to  read  data  files  of  the  form  ENCOM  writes. 
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INFORMATION  ON  OUTPUT  PRODUCTS; 

SETUP:  You  can  see  a  sample  data  file  produced  by  SETUP  at  the  end  of 
Chapter  2«  It  contains  all  the  information  TRBOFN  needs  to  design  and 
calculate  perormance  for  up  to  five  engines  at  the  design  flight  condition 
and  up  to  19  off-design  flight  conditions.  SETUP  also  creates  an  output 
file  used  by  ENCOM  for  input  of  maneuver  conditions  for  sizing  and 
comparison. 

TRBOFN:  The  output  file  from  TRBOFN  includes  design  and  off-design 
information  for  the  engine  simulated  with  information  about  internal 
pressures,  temperatures,  Mach  numbers  and  fuel-to-air  ratios  at  about  50 
stations  inside  the  engine  for  each  operating  condition.  It  also  includes 
thrust  expressed  in  two  ways, thrust  specific  fuel  consumption  and  engine 
efficiencies  at  all  user-input  operating  coonditions.  Thrust  is  expressed 
both  as  a  non-dimensional  quantity  proportional  to  the  physical  size  of 
the  engine  and  a>  thrust  per  total  engine  airflow.  These  expressions  of 
thrust  are  used  in  the  sizing  process  performed  in  ENCOM.  The 
non-dimensional  thrust  expression  is  net  thrust  per  sea  level  static 
pressure  per  flow  area  at  the  first  stage  turbine  vane  choke  point 
(Fq/(PsI  A4i^)).  All  station  areas  in  the  design  process  are 
calculated  as  fractions  or  multiples  of  this  area,  A4]^.  Thus,  the 

non-dimensionalized  thrust  is  directly  proportional  to  engine  size.  In 

other  words,  to  get  a  thrust  of  20,000  lbs  from  an  engine  with  a 
non-dimensional  thrust  of  20.0  would  require  the  product  (Pgi  A41)  to 
be  1,000  and  since  P^l  *  14.69  psi,  A4X  would  be  (1,000/14.69)  sq  in 

or  68.07  sq  in  or  0.4727  sq  ft.  All  the  other  station  flow  areas  are 

calculated  as  fractions  or  multiples  of  A41  so  it  would  then  be  trivial 
to  calculate  any  other  engine  station  flow  area. 

ENCOM:  Two  kinds  of  output  files  are  created  by  ENCOM:  one  kind  contains 
x-y  pairs  for  plotting  with  GRAPH,  the  other  is  a  Mission  Performance 
file.  The  file  for  plotting  is  in  a  simple  form.  The  first  line  in  the 
file  contains  the  number  of  curves  to  be  drawn  on  the  graph,  NLINES.  The 
rest  of  the  file  is  made  up  of  NLINES  groups  of  numbers,  each  in  the  same 
form.  Each  group  consists  of  the  number  of  x-y  pairs  that  are  to  be 
connected  to  make  the  line,  NPTS,  followed  by  that  many  x-y  pairs,  one  on 
each  line  in  the  file.  GRAPH  will  plot  the  numbers  in  any  file  in  this 
form.  The  mission  performance  output  file  contains  performance 
information  on  each  of  up  to  five  engines  for  each  of  the  eleven  legs  of 
the  mission.  It  includes  such  information  as  thrust  required  on  each  leg, 
time  to  complete  the  leg,  distance  travelled,  fuel  used  on  each  leg  as 
fraction  of  take-off  gross  weight  of  the  aircraft  and  wing  loading  used  on 
each  leg  based  on  remaining  fuel  fraction  and  take-off  wing  loading.  The 
bottom  line  output  in  this  file  is  the  minimum  fraction  of  take-off  gross 
weight  that  must  be  fuel  for  the  airplane  to  complete  the  mission  using 
each  engine  under  consideration. 

GRAPH:  Publication  quality  plots  with  user-selected  labels  and  axis 
scaling  can  be  made  with  GRAPH.  Many  examples  of  outputs  from  GRAPH  are 
included  in  this  document,  especially  in  the  Appendix. 
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Chapter  2. 

SETUP:  Interactive  Graphics  Design  Input  Program 
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INTRODUCTIOM; 


SETUP  is  the  interactive  graphics  program  used  in  this  system  to  input  and 
change  engine  design  and  off-design  information  which  is  used  by  TRBOFN  to 
calculate  the  performance  of  the  engine(s)  you  design  at  the  flight 
conditions  you  request* 

SETUP  is  easy  to  run  since  almost  all  the  instruction  you  need  to  run  it 
appears  on  the  screen.  This  chapter  of  the  User's  Guide  is  intended  to 
show  you  each  screen  you  would  see  running  SETUP  on  the  terminal  for  which 
it  was  written,  the  Tektronix  A112.  If  you  don't  have  access  to  a  4112 
terminal,  you  may  have  to  modify  SETUP  before  it  will  work  on  your 
terminals.  SETUP  is  written  using  Tektronix'  Interactive  Graphics  Library 
(IGL).  If  you  do  not  have  access  to  IGL  but  you  do  have  access  to  another 
graphics  library,  you  can  write  a  SETUP  program  of  your  otra  using  your 
graphics  library  to  do  the  things  described  in  this  chapter. 
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When  you  run  SETUP,  Che  firaC  thing  you  aee  is  an  aniaaCad  angina  building 
the  title  screen.  When  Che  aninated  engine  is  finsihed  flying  across  Che 
screen  it  looks  like  Figure  1  below. 


TURBINE  ENGINE 
TECHNOLOQY''PERFORMflNCE 
EUflLUPTIQN  SYSTEM 


r.j.  ariLta  aasaaaai  Lsaoasncni'  tees 


Figure  1.  Title  Screen 

After  a  pause  of  a  few  seconds,  the  screen  goes  blank  and  the  next  screen 

eppears. 

The  first  input  paraaeter  you  get  to  choose  is  Che  number  of  engines  you'd 
like  to  design  on  this  run  through  tlTUF.  later  a  nuaber  froa  I  to  S  and 
press  return  to  get  to  the  next  screen. 


ENTER  1>C  NUTBER  OF  ETBINES  FOR  kMICH  VOJ 
UOULO  LIKE  ID  ENIER  DESIGN  PWMCTERS. 

PLL  DESIGN  DTORmnON  UHL  K  MITIEN  TO  INE  SATC 
DNTR  FILE  TO  RE  RERD  tY  He  PROBNin  TRMFN 


2.  Selecting  Ruaber  ef  Ingiane 


Once  you  enter  the  nuaber  of  enginee  you  went  to  deeign,  the  prograa 
announces  that  you  are  going  to  be  selecting  design  paraaeters  for  engine 
nuaber  1  as  follows: 

ENGINE  DESIGN  PARAMETER  SELECTION 
FOR  ENGINE  NUMBER  1 

IN  -mxe  eccTSON  of  -me  peoeaAM  vou  uxu.  ec  alloueo 
TO  SCLeCT  AND  CHAN6C  ALL  -nC  PARArClCRS  NECCBSARV 

FOR  -me  PROORAn  trsofn  to  OESieri  and  calculate 
oesioii  ANo  -mROTTLeo  FERraeriAncc  for  thc  erMxne 

Figure  3.  Announceaent  Page  for  Engine  Nuaber  1 
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TKBOTM  ia  capable  of  eiiaulating  parforaance  of  both  nixed  and  aeparaca 
flow  ancine  daaigna  ao  the  next  thing  you  get  to  apecify  for  the  engine 
you're  deaigning  ia  which  type  it  ia. 


MIXED  FUOU  ENGINE 


SEPARftTE  FLOU  ENGINE 


ONDOse  HDce  on  scmmie  fum  cmnc  uvaii 

Figure  4.  Hixed  or  Separate  Flow  Oeaign 


The  aext  screen  to  eppeer  is  the  first  design  pareaeter  aenu.  The  next 
two  figures  show  how  highlighting  works.  All  the  pereaeters  on  this  first 
aenu  refer  to  the  coapression  section  of  Che  engine  so  the  coapression 
section  of  Che  engine  drawn  in  the  top  half  of  Che  screen  blinks  on  and 
ff- 


HONGE  UALUE(S) 


COrPRESSIOH  SECTION  PMlArCIERS 

DTRL  PRESSURE  RECOUERV.  OIFFUSER  8.99  9VPA8S  RATIO  1. 

SrFRCSSOR  PRESSURE  RATIO  8.98  FAN  PRESSURE  RATIO  3. 

SSOR  AOlAtATZC  EFFICXCHCV  8.98  FAN  A0XA9ATXC  EFFXdENCV  8. 

NEXT  rcnu 


Figure  5a.  Coapression  Section  Highlighting 


EXT  TEMU 


COrfRESSION  SECTION 

PAR  Arc  TEAS 

rOTRL  PRESSURE  RECOICRV.  OIFFUSER  8.99 

8VPA8S  RATIO 

1.90 

XlraRESSOR  pressure  ratio  9.98 

FAN  PRESSURE  RATIO 

3.99 

MWFRES90R  A02A9ATIC  EFFICtci^  8.98 

FAN  ADIA8ATIC  EPFZCXENCV 

8.98 

NEXT 

ISMJ 

Figure  5b.  Coapression  Section  Highlighting 


To  chance  pereaetere  in  Che  acnu,  aove  Che  croseheir  cursor  onCo  'VTHAWpB 
VAUIBCS)"  end  press  any  alphenusMric  key  or  tCTIHUI  (rigure  5b).  The 
cursor  juaps  inCo  Che  aenu  seccimt  of  Che  screea  and  is  placed  on  "RCXT 
NEHU"  eC  Che  boCCoa  of  che  box  as  shown  in  figure  6. 


XT  «riu 

AN6E  UALUeCS) 


'niNL  Mcssuee  eccoucev,  oxrruscR  a. 
CBiaaceaoa  eactauec  mtzo  s. 


saoe  eoiMarxc  CFricicncv  a. 


Figure  6.  Menu  I:  Coapression  SecClon  Paraaecers 


When  you  press  the  key,  the  value  the  cursor  was  on  disappears  from  the 
menu  and  a  prompt  for  the  parameter  you  chose  appears  at  the  bottom  of  the 
screen.  The  prompt  string  you  see  at  the  bottom  of  the  screen  is  the  muio 
of  the  paraiseter  used  in  the  program  SETUP  and  in  the  engine  simulation 
program,  TRBOFN. 


flail 


Figure  8.  Prompting  for  a  New  Value 


U 


AC  this  point,  encer  Che  new  value  you  want  Che  paraawcer  Co  have  (in  Chia 
example,  2.7)  and  press  RETURN.  The  new  value  will  reappear  in  the  place 
where  the  old  one  was  and  the  cursor  will  reappear  on  “NEXT  HENU".  If  you 
want  to  change  any  other  parameters,  move  Che  cursor  onto  Che  value  you 
want  to  change  and  follow  Che  same  procedure.  If  all  Che  values  on  Che 
menu  are  satisfactory,  simply  press  any  key  while  Che  cursor  is  on  Che 
"NEXT  MENU”  selection  in  the  menu  and  Che  next  menu  will  appear. 


NEXT  rENU 
CHANGE  UALUECS) 


COrS'RESSIOH  SCCTCON  P»«fttCTERS 


TOTW.  PRESSURE  RECOVERY,  DIFFUSER  0.  *9  BYPASS  RATIO 


COMPRESSOR  PRESSURE  RATIO 


S.Se  FAN  PRESSURE  RATIO 


COMPRESSOR  ADIABATIC  EFFICIENCY  e.fe  FAN  ADIABATIC  EFFICIENCY  0.90 


Figure  9.  New  Value  in  Place  of  Old  One 


The  secoad  mnu  iacludea  ell  Che  pareaetera  Che  prograa  TRBOFN  needa  Co 
celculace  Che  effecCa  of  uaiag  an  alCernaCe  fuel  co  run  Che  engine.  The 
paraaeCera  chac  appear  by  defaulc  (ahown  in  Fig  10)  are  for  JP-4.  To 
change  any  or  all  of  chea,  follow  che  aaae  procedure  aa  deacribed  above 
for  changing  che  fan  preaaure  racio.  The  burner  apliccer  place  ia 
highlighced  (i.e.,  blinka  on  and  off)  on  chia  acreen. 


Figure  10.  Menu  2;  AlCernaCe  Fuel  SelecCion 


16 


8ai»*>awwi 


Whea  you  select  either  "NEXT  MENU"  this  screen  appears.  TRBOFN  needs 
several  internal  Mach  numbers  at  design  to  calculate  flow  areas  in  the 
engine  for  use  in  off-design  performance  and  balancing  calculations.  The 
third  menu  allows  you  to  change  four  of  the  necessary  design  Mach 
numbers.  The  station  designator  lines  are  highlighted  on  this  screen. 


Figure  11.  Menu  3:  Internal  Flow  Mach  Numbers 
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the  next  menu  allows  you  to  select  the  control  scheme  to  be  used  to 
determine  the  intermediate  power  setting  for  the  engine  you're  designing 
when  it  runs  at  off-design  flight  conditions.  You  can  choose  to  have 
TRBOFN  consider  either  of  two  temperatures  as  one  of  the  possible  control 
parameters  while  searching  for  off-design  intermediate  power  settings: 
the  turbine  inlet  gas  temperature  or  the  first  stage  vane  metal 
temperature.  If  you  choose  to  have  it  check  gas  temperature,  this  is  the 
menu  where  you  get  to  set  its  limiting  value.  You  can  set  the  maximum 
allowable  metal  temperature  on  a  different  menu. 


■ 

Sil 

IP— ftSB 

tsExr  rtNu 

CHANGE  UALUECS) 


EflGIME  COMTROL  SCHEME 

ftMSUER  THE  OUESTIOM  SEV.OU  TO  DETERMIHE  UHETWR  THE 
ENSIME  UILL  USE  THE  TURBINE  INLET  GAS  TENPERATURE 
OR  FIRST  STAGE  WANE  NETAL  TETP  AS  ONE  OF  THE  POSSIBLE 
CONTROL  PARATCTERS  IM  FINDING  THE  IMTERTtDIATE  POUER 
SETTING  AT  OFF -DESIGN  FLIGHT  COHOITIOHS 

IF  COHTROLLIMG  TO  GAS  TETP,  SET  MAX  T40  GAS:  ZaBV.O  R 


NEXT  MENU 


00  YOU  WANT  THE  ENGINE  TO  BE  CONTROLLED  TO  GAS  TEfPERATURE,  TAB?  lY'HJ: 
Figure  12.  Menu  4:  Control  Parameter  Selection 
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When  you  answer  the  question  "yes",  the  cursor  appears  in  the  menu  block 
on  "NEXT  MENU".  At  this  point  you  can  change  the  value  of  the  gas 
temperature  to  be  used  as  a  control  parameter  from  2220  R  to  whatever  you 
like.  Place  the  cursor  on  "NEXT  MENU"  and  hit  any  key  to  go  on  to  the 
next  menu. 


— 

•* 

NEXT  r^NU 

CHANGE  UALUECS) 

EtIGIME  CONTROL  SCHETt 

ANSUER  THE  OUESTIOM  BELOU  TO  DETERMINE  UHETMER  THE 
ENGINE  LULL  USE  THE  TURBINE  INLET  GBS  TEtTERBTURE 

OR  FIRST  STAGE  UftNE  ftTBL  TEfV»  AS  ONE  OF  THE  POSSIBLE 
CONTROL  PARAMETERS  IN  FINDING  THE  INTERTtDIATE  POLCR 
SETTING  AT  OFF -DESIGN  FLIGHT  CONDITIONS 

ir  rnMTD^  i  Tn  /:ac  tti-d  cf-r  T^a  /:of.  g  f) 

Figure  13.  Menu  5:  Setting  the  Gas  Temperature  for  Control 
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On  the  next  menu  you  get  to  set  the  -adibetic  efficiency  of  the  two 
turbines  in  the  engine.  Change  them  in  the  same  way  as  any  other  menu 
parameter. 


Figure  14.  Menu  6:  Turbine  Adiabatic  Efficiencies 
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iPMMi 


MHnaii 


NeKt,  you  mu«t  answer  another  question.  This  time  you  must  decide  %rhether 
you  want  to  set  any  turbine  cooling  parameters.  For  the  purposes  of 
showing  as  many  of  the  program  screens  as  possible,  let's  answer  the 
question  "yes"  and  see  all  the  possible  turbine  cooling  parameters  you  can 
set. 


OO  YOU  WISH  TO  SET  OMY  TURBinC  COOLING  P«R«NETE»S7  tY-713: 


Figure  15.  Question  on  Setting  Turbine  Cooling  Parameters 


The  iaatrucCioiM  oo  how  to  aoloct  m  high  prosaur*  turbia*  vaa«  coolint 
achaaa  taehoology  laval  ara  aalf-aaplanatory«  load  tha  taat  at  tha  bottoa 
of  aach  of  tha  followiag  acraaaa.  PMCTVH  aCaada  for  farcaat  cora  Maaaflow 
for  Cooliag  tha  ^rbiaa  Vaaaa  of  tha  High  praaaura  turbiaa. 


SPCCIFYINfi  TMC  TECHMOLOGY  LCUEU  OF  n«  HP  TUtaitC  UMC  COOLING  SCHCrC 


THIS  SCLCCTION  IS  MMC  FROn  A  SCRIES  OF  CHARACTER  ISTICAU.Y  SHAPED  CURUCS 
RCPRCSCNTING  TEN  UtLUCS  OF  THE  COOLING  SCHCtC  TECHNOLOGY  PARMCTER.  THE 
CHARACTERISTIC  SHAPE  IS  SIMILAR  TO  A  HtPEROOLlC  TANGENT  CURVE.  THE  COOLING 
SCMCrC  TECHNOLOGY  PARATCTER  IS  THE  COOLING  CFFECTIUCNCSS  OF  1HE  SCHEtC  USING 
ia«  OF  THE  TOTAL  CORE  AIRFLOU  FOR  UANE  COOLING.  SELECT  A  COOLING  SCHETC 
TECHNOLOGY  BY  PICKING  A  POINT  UCRE  PMCTVH'IB*  ANL  THE  VALUE  OF  THE  COOLING 
EFFCCTIUENESS  IS  THE  VALUE  YOU  UANT  TO  USE  FOR  THE  TECHNOLOGY  PARAfCTER. 


NEED  ANY  TORE  INSTWJCTIONST  lY/Hl: 


Figure  16a.  laatructioas  for  Selecting  a  High  Preaaure  Turbine  Vane 
Cooling  Scheme  Technology  Level 
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Figure  16d.  HP  Turbine  Vane  Cooling  Scheme  Technology  Level 


This  neau  includes  all  Che  paraaeCers  necessary  ior  TRBOFN  to  siwilaCe 
cooling  of  Che  high  preaaure  turbine's  vanes  and  rotors.  Mote  Chat  Che 
vane  cooling  technology  parasMter  is  on  this  menu.  That's  the  value  you 
selected  on  Che  last  screen.  The  default  values  for  the  airflow  and 
effectiveness  of  the  rotor  cooling  scheae  are  set  for  no  cooling.  Placing 
seros  in  these  paraaeCers  results  in  the  prograa  crashing.  The  vane  metal 
temperature  allowable  is  only  iaporCant  if  you  chose  aetal  temperature  as 
the  control  parameter  instead  of  gas  temperature. 


E 


dHOMGE  UALUE(S) 


HIGH  PRESSURE  TURBINE  COOLING 
unHE  COOLING  ROTOR  COOLING 

pRNE  COOLING  TECH  PHRRM  0.731  COOLING  8.001 

COOLING  EFFECTILCNESS  0.081 

MACH  NUMBER  OF  MAIN  FLOU  AT  MIXING  PLANE 
hN  FIRST  STAGE  UANES  (4A>  0.800  EXITING  ROTOR  (AD> 


UaNE  rCTAL  TEMP  ALLOUABLE  2880. R 


0.500 
NEXT  MENU 


Figure  18a.  Menu  7:  High  Pressure  Turbine  Cooling  ParaaeCers 


Figures  18a.  and  18b.  show  the  high  pressure  turbine  and  the  two  cooling 
air  ducts  highlighted  in  the  engine  drawing. 


Figure  18b.  High  Pressure  Turbine  and  Cooling  Duct  Highlighting 


The  next  oenu  includea  all  the  parameters  necessary  for  low  pressure 
cooling.  The  only  differences  are  that  the  vane  cooling  scheme  is 
specified  by  a  set  value  of  airflow  and  effectiveness  and  there  is  no 
option  to  set  metal  temperature  allowables  for  the  low  pressure  turbine 
vanes.  Note  the  low  pressure  turbine  and  cooling  ducts  are  highlighted  in 
the  engine  drawing  that  accompanies  this  sienu. 


>T  ^ _ J 

<ICWT  HCHU - 

CHANGE  OALUE(S) 


LOU  PRESSURE  TURBINE  COOLING 

UANE  COOLING 
CORE  RIR  FOR  COOLING 


OOLING  EFFECTIOENESS 


e.BBI 

e.eei 


ROTOR  COOLING 
CORE  AIR  FOR  COOLING 


COOLING  EFFECTIUENESS 


e.BBI 

e.BBI 


MACH  HUMBER  OF  MAIN  FLOU  AT  MIXING  PLANE 
pN  first  stage  vanes  (4R)  B.SBB  EXITING  ROTOR  (4U> 


e.SBB 
NEXT  MENU 


Figure  19.  Menu  8:  Low  Pressure  Turbine  Cooling  Parameters 
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For  mixed  flow  eagine  designs,  you  have  the  option  of  having  augmented 
performance  calculated.  Separate  flow  designs  can  not  be  augmented  with 
TRBOFN.  If  you  do  want  augmented  performance  calculated  for  a  mixed  flow 
design  then  as  soon  as  you  answer  this  question  "yes”  you  will  see  the 
next  question  which  appears  with  Figure  21. 


NEXT  MENU 
CHANGE  UALUE(S) 


00  YOU  WISH  TO  aoCrCMTED  PERFORMBMCE  CaLCULSTtO?  tY^N3i 


Figure  20.  Augmented  Performance  Option 


The  nuoiber  of  augmented  performance  points  refers  to  throttling  the 
afterburner.  Entering  one  point  tells  the  program  to  simulate  only  full 
afterburner;  any  larger  number  tells  it  to  first  turn  on  the  afterburner 
full  and  then  throttle  it  back  to  a  near  off  condition  in  that  many  steps. 


ErntR  THE  nureER  or  AUCnEnTED  points  for  wich  you  uduld  like 
PERFORMBNCE  CftLCULRTED  «T  THE  DESIGN 
FLIGHT  CONDITION:  Ce-91 


Figure  21.  Number  of  Augmented  Points 


The  next  oeau  allows  you  to  specify  how  much  of  the  total  engine  airflow 
is  to  be  used  for  nozzle  liner  and  flap  cooling.  TRBOFN  considers  these 
cooling  ducts  to  be  open  for  all  operating  points  for  any  augmented  engine 
with  nozzle  cooling.  The  default  values  are  for  no  cooling.  Once  again, 
entering  zeros  for  these  values  results  in  division  by  zero  in  TRBOFN. 
The  cooling  effectivenesses  are  for  use  in  calculating  the  resultant 
material  temperatures  of  the  nozzle  liner  and  flap. 


Figure  22.  Menu  9:  Nozzle  Cooling  Parameters 
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The  next  screen  is  just  an  announcement  page  to  tell  you  that  next  you  get 
Co  select  Che  flight  conditions  where  you  want  performance  calculated. 


SELECTION  OF  DESIGN  POINT 

AND 

POINTS  FOR  OFF-DESIGN  PERFORMANCE 
FOR  ENGINE  NUMBER  1 


Figure  23.  Flight  Condition  Selection  Announcement 
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In  specifying  flight  conditions  for  the  engine  performance  program  to  use, 
you  can  specify  directly  the  number  of  throttled  power  settings  to 
calculate  or  you  can  tell  the  program  to  use  the  flight  conditions  as 
maneuver  requirements  in  which  case  the  program  will  automatically  set  the 
number  of  throttled  power  settings  at  which  to  calculate  performance  at 
the  maximum  value  of  19.  The  maneuver  requirement  option  is  useful  if  you 
plan  to  compare  several  engines  with  the  engine  sizing  and  comparison 
program,  ENCOM.  ENCOM  compares  engines  by  requiring  them  to  propel  an 
airplane  through  a  series  of  maneuvers  which  are  specified  as  sustained 
accelerations  at  specific  Mach-altitude  combinations.  If  you  specify  the 
maneuver  requirements  you  intend  to  use  in  ENCOM  now,  you  can  be  sure  your 
engine  performance  data  files  will  include  the  performance  information 
ENCOM  will  need  to  do  its  comparison.  Alternatively,  if  you  don't  want  to 
use  ENCOM,  you  can  specify  the  number  of  throttled  power  settings  directly 
for  each  flight  condition  you  choose.  This  will  make  TRBOFN  run  faster 
since  it  will  have  fewer  performance  points  to  find  and  calculate. 


TO  SPECIFV  OFF-OESIGN  FLIOHT  COHOITIONS  IKERE  VOU  UJULD  LIKE  PERFDRfttNCE  CflLCULflTEO 

YOU  HflUE  the  options  OF  ENTERIMB  EITHER: 

1.  mnEULAR  REOUIREIWTS:  (MRCH,  HIT,  SUSTftIMED  C  CAPABILITY) 

DR 

P.  nUCH-ALTiniDE  PAIRS  UITH  flUr*ER  OF  THROTUEO  POk£R  SETTINGS  AT  EACH 

EfTTER  1  roe  MAMElACR  REOUIREICNTS;  2  FOR  PnCH-RLT  PAIRS: 

Figure  24.  Maneuver  Requirements  or  Mach-Altitude  Pairs  and  Number  of 
Throttled  Power  Settings  at  Each 
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In  either  case,  the  next  screen  you  see  after  answering  the  question  is  a 
plot  of  Altitude  vs  Mach  number.  You  can  change  the  limits  on  the  axes  to 
suit  your  needs  very  simply  by  answering  the  questions  as  in  the  sequence 
depicted  in  Figures  25a  through  c.  Hitting  "&£TURN"  trtien  the  program  asks 
whether  you'd  like  to  change  the  limits  on  one  of  the  axes  is  the  sasu  as 
answering  "no." 


MODEL  COMSTRaitiTC: 

itLTnuoc  e  TO  iBs.aae  ft 
tKCH  HO.  e  TO  8.S 


Figure  25a.  Changing  Limits  of  Altitude  vs  Mach  Number  Space 


34 


n.i(jHT  rmcH  nutser 


UOUt-O  rtXJ  LIRE  TD  CHANGE  LIMITS  OF  MACH'>Y 
EnrER  MACH  NO.  MIN  THEN  MAX: 0,2 

Figure  25b.  Changing  Limits  of  Altitade  vs  Mach  Number  Space 


UOULO  you  LIKE  TO  CHANGE  LIMITS  OF 
UOULD  you  LIKE  TO  CHANGE  LIMITS  OF  MACH? 


NUrOER 


Figure  25c.  Changing  Limits  of  Altitude  vs  Mach  Number  Space 


When  you're  saClsfied  with  the  limits  on  the  axes  you  get  to  use  the 
crosshair  cursor  to  select  the  design  flight  condition  for  the  engine. 
Use  the  thumb  wheels  on  the  4112  terminal  to  move  the  cursor  to  the 
Hach-altitude  you  want  and  press  any  key.  The  Hach-altitude  pair  you 
picked  appears  next  to  the  marked  point  you  chose.  In  addition,  a  line  of 
constant  fan  inlet  total  temperature  (■'^20^  drawn  through  the  point 
you  picked.  This  is  an  aid  in  selecting  high  and  low  T2Q  points  for 
off-design  operation  and  for  doing  a  study  at  constant  T^g.  If  the 
values  are  not  what  you  thought  they  would  be,  you  can  type  in  the  values 
you  intended  by  answering  the  question  in  the  dialog  area. 


noKL  CDNsrnaiMTS! 

M.TnUOC  0  TO  IBS, 
MBCH  NO.  a  TO  2.5 


FT 


use  CIKSQR  TO  PICK  OCSIGN  POINT 


t- 

U. 

§ 

P 


n.i«T  mcM  NuniCR 

UDULO  VtJU  LIKC  TO  CHWIGE  LIMITS  OT  MOCH-> 

IS  TMaT  OMT  YOU  INTtNOCO'*  tV^NI: 


Figure  26.  The  Design  Point  and  Design  Fan  Inlet  Total  Temperature 
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When  you  have  selected  the  design  point  you  intended,  you  get  to  enter 
either  the  number  of  throttled  power  settings  you  want,  or  the  sustained 
acceleration  you  want  the  aircraft  to  hold  at  this  flight  condition, 
depending  on  whether  you  chose  to  enter  maneuver  requirements  or  not. 


EMTES  T>c  MUMJER  OF  TOROTTV-EO  POUER  SETTINGS  FOR  (JHICH  yOU  UOULO  LIKE 
PERFORmoKE  CflLCULRTEO  AT  THIS  FLIGHT  CONDITION:  10-19]: 


Figure  27.  Number  of  Throttled  Power  Settings 


If  you  want  to  fly  the  engine  through  the  mission  in  ENCOH,  you  need  to 
have  TRBOFN  calculate  performance  for  the  engine  at  the  flight  conditions 
used  in  the  mission  profile.  SETUP  will  build  the  input  file  to  make 
TRBOFN  run  the  engine  at  all  the  necessary  flight  conditions  and  it  will 
also  make  a  mission  input  file  to  tell  ENCOM  the  maneuver  requirements  you 
want  the  engine  to  meet  if  you  answer  "yes"  to  the  question  in  Figure  28. 


Figure  28.  To  Build  Mission  Input  File  for  ENCOM 
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Either  way,  you  get  the  option  to  select  more  flight  conditions  if  you 
want  to. 


e.a  B.s  i.a  i.s  2. a 

FLIGHT  mCH  HUrBCR 


DO  YOU  UflllT  TO  FLY  THIS  ENGINE  THROUGH  THE  MISSION  IN  ENCOM?  CY'NJiN 
UftiT  TO  RICK  OTHER  FLIGHT  CONDITIONS  FOR  THIS  ENGINE?  CY/'NJ: 

Figure  29.  More  Flight  Conditions 


Selecting  these  off-design  flight  conditions  is  just  like  selecting  the 
design  point  in  that  you  move  the  cursor  to  the  desired  point,  hit  a  key 
to  encode  it,  and  the  program  displays  the  coordinates  of  the  point  you 
picked  with  a  line  of  constant  T2Q  through  it  and  asks  you  if  that  was 
what  you  intended. 


PICK  UP  TO  19  MORE  FLIGHT  COHDITIOtlS  (OUT  OF  BOUMDS  TO  STOP) 


FLIGHT  MUCH  HUrSER 

URHT  TO  PICK  OTHER  FLIGHT  CONOITIONS  FOR  THIS  EtlGIIIE?  CY^HJiY 
IS  THRT  UWT  YOU  IMTEMOEDT  tY^HJ:Y 


Figure  30.  Picking  an  Off-Design  Flight  Condition 


When  you  are  satisfied  with  the  flight  condition  on  the  display,  you  get 
to  enter  the  number  of  throttled  qugmentor  points  you  want  calculated  if 
your  engine  is  a  mixed  flow,  augmented  design. 


PICK  UP  TO  19  ItKC  FLIGHT  COHDITIONS  (OUT  OF  BOUflOS  TO  STOP) 


ENTER  THE  NurCCR  OF  AUGrCNTEO  POINTS  FOR  LMICH  YOU  UOULD  LIKE  PERFORMANCE 
CRLCULATED  AT  THIS  FLIGHT  CONDITION: 
ce'9]  N0TE:C1  FOR  MAX  AS  ONLY!: 


Figure  31.  Number  of  Throttled  Augmentor  Points 


When  you  have  chosen  all  the  off-deaign  points  you  want  (up  to  19),  pick  a 
point  outside  the  limits  of  the  Mach-altitude  axes  and  the  program  will 
erase  the  graphics  screen  and  ask  you  for  the  name  you  want  the  data  file 
it  makes  to  have.  This  is  the  file  that  TRBOFN  reads  to  design  and 
calculate  performance  for  the  engine  you  just  described. 


The  output  file  from  SETUP  is  the  input  file  for  TRBOFM.  Below  is  an 
example : 
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T 


U.94U000U 

1 .000000 

3.0O0P0O 

O.VOOOoOo 

b .000000 

O.9OUOO0O 

1 35.0y0o 

lHt>no.uO 

O.SOOnOOo 

51b.bPil0 

6.77b233tfr-02 

0.4000000 
7*20. yOo 
O.VOoOOOo 

o.'ioCoOof 

k2u0.0o0 

C.OuOoOyO 

0.1000000 

T 

0. 1 bOyO JO 

0.5002206 

O.HoOOPoP 

l,0P0O000F-O3 

i.ouooooct-os 

0.5000000 

1  ,ij0  jOOOuF-03 
O.5000000 

1  . 000000011-03 

O.hOOOuOq 

1 .OyOyOOOt-03 

l.OOOOOOuF-03 

T  !(.v9oOv«7r- 

c 

(■'5  o.oo9y®97t- 

o5  0.5(i0<io0C 

O.OoOyOyOt^OO 

J 

1  .'oil?/. 

1007l .OK 

10 

9 

0.»97P021 

1C 

9 

o.uOoOoOcrtoo 

o.roOoOoOU+OP 

10 

C 

1 .703730 

11920. vt 

10 

9 

0.  /95S34i 

4'.0«3.o0 

lO 
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The  variables  that 

these  numbers 

represent  are: 

MEMG 

MIXED 

PQPOIF 

FPR 

FEA 

CPR 

CEA 

BR 

NWF 

FLHV 

CPF 

TOF 

FARSTO 

XM20 

XM26 

XM33 

XM50 

T40MVA 

T40G 

CONGAS 

THEA 

TLEA 

CVTECH 

XM4A 

PMCTRH 

ECTRH 

XH4D 

PMCTVL 

ECTVL 

XM4R 

PMCTRL 

ECTRL 

XH4U 

AUGMNT 

PMCNL 

PMCNF 

ECNL 

ECNF 

NPAIRS 

XH(1) 

ALT(l) 

NPWR( 1 ) 

NAUGPT(l) 

XM(2) 

ALT(2) 

NPWR(2) 

NAUGPT(2) 

• 

• 

a 

• 

• 

• 

• 

• 

• 

• 

• 

• 

•  •  •  • 

XM(NPAIRS)  ALT(IIPAIRS)  NPWR(NPAIRS)  liAUGPT(NPAIRS) 


42 


The  paremeters  that  appear  in  the  data  file  are  in  the  sane  order  as  you 
entered  them  while  running  SETUP.  The  following  list  defines  all  the 
variable  names  listed  on  the  last  page. 

NENG  Number  of  ENGines  for  which  there  is  design  and  off-design 
information  in  this  file 

MIXED  Logical  flag  to  denote  mixed  (T)  or  separate  (F)  flow  design 

PQPOIF  Fraction  of  total  pressure  lost  across  diffuser  due  to  friction 

FPR  Fan  Pressure  Ratio 

FEA  Fan  Efficiency,  Adiabatic 

CPR  Compressor  Pressure  Ratio 

CEA  Compressor  Efficiency,  Adiabatic 

BR  Bypass  Ratio  (Note:  TRBOFN  calculates  BR  for  mixed  flow 

designs  based  on  total  pressure  match  between 
the  core  and  bypass  streams  at  the  mixing 
plane  so  the  value  you  enter  for  BR  for  a 
mixed  flow  design  is  ignored  by  TRBOFN) 

MWF  Molecular  Weight  of  the  Fuel 

FLHV  Fuel  Lower  Heating  Value  (Btu/lb) 

CPF  Specific  Heat  of  Fuel 

TOF  Entering  Temperature  of  Fuel 

FARSTO  Stoichiometric  Fuel-to-Air  Ratio 

T40MVA  Maximum  allowable  metal  temperature  of  high  pressure  turbine's 
first  stage  vanes 

T40G  Maximum  allowable  gas  temperature  entering  high  pressure  turbine's 
first  stage  vanes 

CONGAS  Logical  flag:  (T)  for  gas  temperature  T40G  as  control  limit 

(F)  for  metal  temperature  T40MVA  as  control  limit 
THEA  High  pressure  turbine's  adiabatic  efficiency 
TLEA  Low  pressure  turbine's  adiabatic  efficiency 
CVTECH  High  pressure  turbine  vane  cooling  scheme  technology  level 
parameter 

XM4A  Mach  number  of  main  stream  at  mixing  plane  for  high  pressure 
turbine  vane  cooling  air  input 

PMCTRH  Percent  core  Massflow  used  for  Cooling  the  Rotors  of  the  High 
pressure  turbine 

ECTRH  Effectiveness  of  the  Cooling  of  the  Turbine's  Rotors,  High 
pressure 

XM4D  Mach  number  of  main  stream  at  mixing  plane  for  high  pressure 
turbine  rotor  cooling  air  input 

PHCTVL  Percent  core  Massflow  for  Cooling  of  the  Turbine's  Vanes,  Low 
pressure 

ECTVL  Effectiveness  of  the  Cooling  of  the  Turbine's  Vanes,  Low 
pressure 

XM4R  Mach  number  of  main  stream  at  mixing  plane  for  low  pressure 
turbine  vane  cooling  air  input 

PMCTRL  Percent  core  Massflow  used  for  Cooling  the  Rotors  of  the  Low 
pressure  turbine 

ECTRL  Effectiveness  of  the  Cooling  of  the  Turbine's  Rotor's,  Low 
pressure 
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XM4U 


AUGMNT 

PMCNL 

PriCNF 

ECNL 

ECNF 

NPAIRS 

XM 

ALT 

NPWR 

NAUGPT 


Mach  number  of  main  stream  at  the  mixing  plane  for  low  pressure 
turbine  rotor  cooling  air  input 

Logical  flag:  (T)  for  augmented  engine,  (F)  for  non-augmented 
engine 

Percent  total  engine  Maasflow  used  for  Cooling  the  Nozzle  Liner 

Percent  total  engine  Maasflow  used  for  Cooling  the  Nozzle  Flaps 

Effectiveness  of  the  Cooling  of  the  Nozzle  Liner 

Effectiveness  of  the  Cooling  of  the  Nozzle  Liner 

Number  of  Mach-altitude  Pairs  (flight  conditions)  at  which  the 

user  wants  performance  calculated 

Array  of  Mach  numbers  of  flight  conditions 

Array  of  altitudes  of  flight  conditions 

Array  of  number  of  throttled  power  settings  to  be  calculated  at 
each  flight  condition 

Array  of  number  of  throttled  augmentor  settings  to  be  calculated 
at  each  flight  condition 
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TRBOFN: 


Chapter  3. 

The  Ideal  Gas  Engine  Performance  Simulation  Program 
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INTRODUCTION: 


1.  Background:  TRBOFN  was  developed  in  a  project  at  the  Frank  J.  Seiler 

Research  Laboratory  at  the  Air  Force  Academy.  One  of  the  goals  of  the  project 
was  to  develop  a  method  of  assessing  the  relative  benefits  of  various 
improvements  in  the  technology  level  of  the  components  of  a  turbofan  engine 
which  was  fulfilled  by  the  system  documented  in  this  guide.  The  original 
engine  performance  simulation  model  used  in  the  project  was  based  on  perfect 
gas  assumptions  including  constant  specific  heats.  This  proved  to  be  too 
simple  a  model  for  use  in  this  facet  of  the  work.  A  much  more  realistic  ideal 
gas  model  in  which  the  specific  heat  of  the  gas  is  a  function  of  its 
constituents  and  its  temperature,  consistent  with  the  model  used  by  Keeenan, 
Chao,  and  Kaye  in  generating  Gas  Tables  is  used  in  TRBOFN,  the  engine 

simulation  program  described  here. 

2.  Capabilities: 

a.  Design  Parameter  Selection:  Performance  of  almost  any  reasonable 
cycle  turbofan  engine  can  be  simulated  with  TRBOFN.  Cycle  parameters,  a 
design  point  (Mach~altitude  pair),  and  up  to  nineteen  off-design  flight 
conditions  can  be  set  by  the  user  along  with  several  other  parameters.  Some 
of  the  other  user  selectable  items  include: 

(1)  either  mixed  flow  or  separate  flow  engine  design, 

(2)  whetner  the  engine  is  to  be  controlled  to  gas  temperature 
exiting  the  burner  or  to  metal  temperature  in  the  first  stage  turbine  vanes 
(Limits  for  these  control  temperatures  are  set  by  the  user.), 

(3)  use  of  an  alternate  fuel  (The  default  fuel  is  JP4.), 

(4)  adiabatic  efficiencies  for  the  fan,  compressor  and  each  turbine, 

(5)  fractions  of  core  airflow  to  be  used  as  cooling  flows  in  the 
high  and  low  pressure  turbines'  vanes  and  rotors, 

(6)  cooling  effectiveness  and  mainstream  Mach  numbers  for  each  of 
the  four  stations  where  cooling  air  is  mixed  with  the  main  flow, 

(7)  whether  augmented  performance  is  to  be  calculated, 

(8)  nozzle  liner  and  flap  cooling  flow  fractions  and  cooling 
effectivene..£e8  for  mixed  flow,  augmented  engines. 
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b.  Off-Design  Performance  Envelope:  Up  to  nineteen  off-design  Mach 
number-altitude  combinations  can  be  selected  by  the  user  as  points  for  which 
engine  performance  information  is  to  be  calculated.  Performance  at  the 
intermediate  power  setting  and  at  each  of  up  to  nineteen  throttled  power 
settings  is  calculated  at  the  design  flight  condition  and  at  each  of  the 
off-design  points  selected. 

3.  Constraints  and  Limitations:  Certain  combinations  of  fan  pressure 
ratio,  compressor  pressure  ratio  and  bypass  ratio  may  result  in  the  design  of 
an  engine  which  is  incapable  of  operating  at  certain  flight  conditions.  No 
warning  of  this  is  given  by  the  program  until  it  fails.  Selection  of 
reasonable  engine  cycles  and  operating  envelopes  should  preclude  this  but  it 
is  possible  to  request  the  program  to  simulate  performance  of  an  engine  that 
will  not  operate  at  the  points  you  requested.  In  addition  to  this  mode  of 
failure,  there  are  certain  limitations  to  which  this  program  is  subject  due  to 
the  assumptions  built  into  the  models.  In  general,  the  constraints  and  limits 
of  the  simulation  program  are  as  follows: 

a.  Gas  Temperature:  No  account  for  dissociation  has  been  made  in 

the  ideal  gas  model.  Thus,  selection  of  gas  temperatures  or  metal  control 

2 

temperatures  that  result  in  gas  temperatures  above  about  5000  R  in  the 
burner  will  result  in  an  inaccurate  representation  of  the  gas  properties  at 
stations  where  the  temperature  is  at  or  near  this  level. 

b.  Flight  Mach  Number:  For  simulation  of  engine  performance  at  Mach 
numbers  above  about  2.5  this  may  not  be  a  reasonable  model  to  use  because  it 
does  not  include  an  option  to  account  for  overboard  bleed  air  and  all  engine 
walls  are  assumed  to  be  adiabatic. 


c.  Flight  Altitude:  The  limits  of  the  atmosphere  model  in  TRBOFN 
are  the  same  as  the  limits  of  the  U.S.  Standard  Atmosphere,  1962.  Standard 
lapse  rates  for  temperature  and  pressure  are  used  for  calculation  of  static 
properties  at  given  altitudes  and  atmospheric  discontinuities  are  consistent 
with  U.S.  Standard  Atmosphere  tables.  Performance  at  any  altitude  contained 
in  the  U.S.  Standard  Atmosphere  Tables  can  be  simulated. 


THE  IDEAL  GAS  MODEL 
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THE  IDEAL  GAS  MODEL;  As  mentioned  in  ‘the  Introduction,  calculation  of  the 
thermodynamic  properties  of  the  constituents  of  air  and  the  products  of 
combustion  are  performed  in  the  manner  outlined  by  Keenan,  Chao,  and  Kaye  in 
reference  1.  Several  constants  for  the  constituent  gases  were  taken  from 
references  3,  4,  and  5.  The  approach  is  based  on  the  assumption  that  the 
mechanics  of  the  gas  molecules'  energy  storage  can  be  described  statistically 
based  on  the  make-up  of  the  gas  and  its  temperature.  The  following  sections 
describe  this  approach  very  briefly.  A  detailed  description  of  all  the 
calculations  that  result  in  the  ideal  gas  model  is  contained  in  Keenan,  Chao, 
and  Kaye's  book  and  each  calculation  is  also  documented  in  the  CAS  model 
subroutines. 


1.  The  Make-up  of  the  Gas  in  the  Engine:  Air  is  taken  to  be  a  mixture  of 
the  following  gases  in  the  stated  proportions: 

(Nitrogen)  78.032  by  volume 

O2  (Oxygen)  20.99%  by  volume 

Ar  (Argon)  0.98%  by  volume 

The  proportions  of  the  products  of  combustion  CO^  (Carbon  Dioxide)  and 
H2O  (Water  Vapor)  at  any  given  station  in  the  engine  depend  on  the 

fuel-to-air  ratio  at  the  station  of  interest  which  depends  on  the  amounts  of 

dilution  and  cooling  air  mixed  into  the  flow  downstream  of  the  stoichiometric 

combustion  section  of  the  burner. 


2.  The  Gibbs-Dalton  Relation:  The  Gibbs-Dalton  Relation  simply  states 
that  the  properties  of  a  mixture  of  several  different  gases  can  be  calculated 
as  the  sum  of  the  mole  fraction  weighted  contributions  of  the  constituents. 
This  relation  is  used  in  the  calculation  of  the  thermodynamic  properties  of 
various  air/products-of-combustion  mixtures.  The  molecular  weights  of  the 
constituents  of  air  and  the  products  of  combustion  follow: 

Molecular  Weight  (gm/mole) 


Ar  (Argon)  39.9480 

Nj  (Nitrogen)  28.0134 

0^  (Oxygen)  31.9988 

CO2  (Carbon  Dioxide)  44.01 

H2O  (Water  Vapor)  18.01534 
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3.  Calorically  Imperfect  Gaa:  A  calorically  imperfect  gas  obeys  the 
equation  of  state 

P  V  RT  where  P  ~  static  pressure 

s  s  s 

V  *  volume 

R  "  Universal  Gas  Ck>nstant 
~  static  temperature 

which  makes  it  a  thermally  perfect  gas.  A  calorically  imperfect  gas  also  has 
specific  heats  (at  constant  pressure,  c^;  and  at  constant  volume,  c^) 
which  are  functions  of  temperature  only  with 


c  -  c 
P  V 


R. 


4.  Gases  Made  up  of  Monatomic  Molecules:  Argon  is  the  only  constituent 
gas  made  up  of  monatomic  molecules  whose  properties  are  accounted  for  in  the 
gas  modeling  subroutine  MONATOM.  The  functions  relating  c  and  c  to 

p  V 

temperature  for  gases  made  up  of  monatomic  molecules  are  derived  from  the 

classical  theory  of  equipartition  of  energy.  This  theory  states  that  equal 

amounts  of  energy  are  stored  in  each  of  the  mechanical  processes  available  to 

the  molecules  to  store  energy,  provided  a  large  enough  number  of  molecules  is 

considered.  Each  mechanical  process  is  called  a  degree  of  freedom.  Monatomic 

gas  molecules  have  three  degrees  of  freedom,  all  translational.  The 

equipartition  of  energy  theory  leads  to  c^  and  c^  being  independent  of 

temperature  which  implies  that  caloric  imperfections  for  monatomic  gases  are 

small.  This  agrees  well  with  experimental  determination  of  the  specific  heats 

2 

for  gases  made  up  of  monatomic  molecules. 


5.  Gases  Made  up  of  Diatomic  Molecules:  The  properties  of  gases  made  up 
of  diatomic  molecules  (e.g.,  0^  and  N^)  are  based  on  a  statistical 
thermodynamic  model  which  uses  energy  storage  distributions  based  on 
temperature.  Diatomic  molecules  have  six  degrees  of  freedom.  Consider  the 
dumbell  shaped  molecular  structure  in  Figure  1. 
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Figure  1.  Diatomic  Molecule's  Six  Degrees  of  Freedom 

Its  six  degrees  of  freedom  include  three  translational  modes  (along  the  x,y 
and  z  axes),  two  rotational  modes  (about  the  x  and  z  axes;  the  moment  of 
inertia  of  a  molecule  with  this  structure  about  its  longitudinal  axis  is 
negligible  compared  to  its  moment  of  inertia  about  either  of  the  other  two 
normal  axes)  and  one  vibrational  mode.  An  anharmonicity  correction  term  i's 
included  to  correct  for  errors  introduced  by  assuming  the  molecule  acts  as  a 
rigid  rotator  and  harmonic  oscillator.  The  distribution  of  energy  stored  by 
each  of  these  six  degrees  of  freedom  depends  on  the  temperature  at  which 
thermodynamic  properties  are  to  be  calculated.  At  low  temperatures,  the 
contributions  of  rotation  and  vibration  are  negligibly  small  compared  to  the 
contributions  of  translational  modes.  At  moderate  temperatures,  the 
contribution  of  vibration  to  energy  storage  is  atill  amall  while  rotational 
modes  are  dominant  and  translational  modes  are  fully  active.  The 
characteristic  temperatures  at  which  new  modes  become  active  have  been  found 
empirically  and  documented  in  reference  5.  The  model  for  oxygen  also  accounts 
for  electronic  excitation  above  a  characteristic  temperature,  also  taken  from 
reference  5. 
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6.  Gases  Made  up  of  Triatomic  Moelculea;  The  most  complex  gas  molecules 
accounted  for  in  this  simulation  are  triatomic  H2O  and  The 
statistical  quantum  mechanical  approach  to  calculating  thermodynamic 
properties  of  triatomic  gases  was  viewed  as  too  complicated  and  poorly 
documented  to  pursue.  We  opted  instead  to  use  the  equations  to  account  for 
the  contributions  of  translational,  rotational  and  vibrational  energies  and  to 
add  to  them  an  appropriate  contribution  due  to  electronic  excitation  so  that 
the  properties  calculated  would  match  the  tabulated  data  in  the  Gas  Tables 
book.  In  using  these  equations  to  get  the  best  possible  curve  fits  we  used 
the  documented  values  for  characteristic  temperatures  for  rotation  and 
vibration  for  each  molecule.  Using  this  much  of  the  model  we  could  match  the 
portion  of  each  curve  below  the  higher  of  these  two  temperatures  very  well. 
To  model  the  rest  of  the  curve,  however,  we  had  to  simplify  the  accounting  for 
electronic  excitation.  The  simplification  involved  treating  all  electronic 
excitation  modes  as  a  group  and  iterating  to  find  an  artificial  characteristic 
temperature  to  give  the  most  accurate  curve  fits  for  the  tabulated 
properties.  Figure  2  shows  representations  of  the  structure  of  both  water 
vapor  and  carbon  dioxide  molecules. 


Figure  2.  Water  Vapor  and  Carbon  Dioxide  Molecular  Structures 

Note  that  the  carbon  dioxide  molecule  is  linear  and  therefore  has  three 
translational  and  two  rotational  degrees  of  freedom.  Water  vapor,  on  the 
other  hand,  has  three  translational  and  an  additional  third  rotational  degree 
of  freedom  due  to  its  non-linear  structure.  Figure  3  shows  the  close 
agreement  between  the  tabulated  data  from  Keenan,  Chao  end  Kaye  end  the  curve 
fits  we  generated. 
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Figure  3.  Agreenent  Between  Tabulated  Data  and  Curve  Fite 
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EMGIME  COMPONENT  MODELS:  TRBOFN's  prograoi  structure  is  modular.  That  is, 
each  engine  component  is  modeled  in  a  separate  subroutine.  The  purpose  of 
this  section  is  to  present  information  regarding  the  development  of  the 
component  model  subroutines. 

1.  Overview  of  the  Engine  Model;  Figures  Aa  and  b  are  schematic  drawings 
of  the  two  types  of  engines  whose  performance  can  be  simulated  by  TRBOFM. 
Figure  Aa  depicts  a  mixed  flow  engine;  Figure  Ab,  a  separate  flow  engine. 
Following  the  drawings  is  a  list  describing  the  station  number  designations. 


Figure  Aa.  Schematic  of  a  Mixed  Flow  Engine 


Figure  Ab.  Schematic  of  a  Separate  Flow  Engine 
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STATION  NUMBER  DESIGNATIONS; 


20-Seriea  CompreBBion  System 

20  Lp  CompreBBor  Entrance 

23  Lp  Compressor  Exit 

24  Lp  Compressor  Exit  after  Bypass  Diverted 

26  Hp  Compressor  Entrance 

29  Hp  Compressor  Exit 

30-Series  Heat  Addition  System 

30  Entry  after  LPT  Cooling  Airs  Diverted 

31  Entry  after  HPT  Cooling  Airs  Diverted 

33  Primary  Burner  Entry  after  Dilution  Air  Diverted 
and  Primary  Diffused 

34  Primary  Burner  w/Fuel  Mixed 

36  Primary  Burner  Burned-stoich 

3D  Dilution  Air  Entry 

38  Burner  Exit  -  Mixed  Primary  and  Dilution 

40-Seriea  Turbine  System 

40  Hp  Vane  Entrance 

4A  Mainstream  at  Vane  Cooling  Air  Entry 

4B  Cooling  Flow  at  Vane  Cooling  Air  Entrance 

4C  Mixed-Mainstream  and  Cooling 

41  Choke  Point 

42  Hp  Rotor  Entrance 

4D  Mainstream-Blade  Cooling  Entrance 

4E  Cooling  Air-Blade  Cooling  Entrance 

4F  Mixed-Mainstream  and  Cooling 

46  Lp  Vane  Entrance 

4R  Mainstream-Vane  Cooling  Entrance 

4S  Cooling  Air-Vane  Cooling  Entrance 

4T  Mixed-Mainstream  and  Cooling 

47  Choke  Point 

48  LP  Rotor  Entrance 

4U  Mainstream-Blade  Cooling  Entrance 

4V  Cooling  Air-Blade  Cooling  Entrance 

4W  Mixed-Mainstream  and  Cooling 

4X  After  Fan  ID  Lp  Work  Extracted 

4Y  After  Other  Engine  Related  Lp  Work  Extracted 

4Z  After  Fan  OD  Lp  Work  Extracted 

49  Lp  Turbine  Exit-After  Customer  Lp  Work  Extracted 


SQ-Seriea  Mixer  System 


50  Core  Entrance 

15  Bypass  Entrance 

59  Mixer  Exit 


60-Seriea  Augmentor  System 
60  Entrance 

63  Primary-Dilution  Air  Diverted 

64  Primary-Mixed  w/Fuel 

66  Primary-Burned 

6D  Dilution  Air  Return 

67  Mixed  Mainstream  also  Liner  Cooling  Air  Return 

6E  Coaling  Air  -  Liner  Cooling  Air  Return 

68  Mixed 

69  Exit 

70,  80,  90-Serie8  Nozzle  System 

70  Nozzle  Entrance 

80  Throat 

8A  Mainstream  at  Liner  Cooling  Air 

8B  Cooling  Air  -  Liner  Cooling  Air 

8C  Mixed  Mainstream 

90  Exit  Plane 

lO-Series  Bypass  System 

15  Bypass  Return  to  Mixer  if  Mixed  Flow  Engine 

18  Bypass  Nozzle  Throat  if  Separate  Flow  Engine 

19  Bypass  Nozzle  Exit  if  Separate  Flow  Engine 
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Several  smaller  seccion  schematics  for  engine  components  are  included  in 
the  following  sections  to  illustrate  the  modeling  assumptions  for  each 
component . 


1.  The  Diffuser:  The  total  pressure  recovery  through  the  diffuser 
accounts  for  two  loss  mechanisms: 

a.  Friction:  The  total  pressure  loss  due  to  frictional  effects 
is  accounted  for  as  a  user-selected  fractional  pressure  drop  assumed  to  be 
constant  with  Mach  number  and  altitude. 

b.  Shock  Losses:  The  total  pressure  loss  due  to  shocks  is 
accounted  for  in  accordance  with  MIL  5008-C: 

P3  *  (Pj)  (1  “  ( .075)(M-1)^*^^)  where  P^  “  pressure  after  loss 

P^  •  pressure  before  loss 
M  "•  flight  Mach  number 

for  flight  Mach  numbers  greater  than  1.0. 

These  calculations  are  performed  in  the  subroutine  RAMREC. 


0  8  3 


Figure  5.  Diffuser  Model 
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Calculation  of  other  total  thermodynamic  properties  after  diffusion  based 
on  free  stream  statics  and  flight  Mach  Number  is  performed  in  the 
subroutine  RAM. 


2.  The  Fan  and  Compressor:  The  compression  processes  performed  by 
the  fan  and  compressor  components  are  assumed  to  be  adiabatic.  Thus,  each 
component  has  an  adiabatic  efficiency  associated  with  it  which,  for  a 
compression  process,  is  defined  as  the  ratio  of  the  enthalpy  required  for 
an  isentropic  compression  process  (i.e.,  at  constant  entropy,  "s")  to  that 
required  to  perform  the  same  compression  adiabatically.  These  processes 
can  be  described  by  paths  on  a  plot  of  enthalpy  versus  entropy,  an  H-s 
diagram.  Paths  1  and  2  on  Figure  4  are  lines  of  constant  pressure  where 
Pj  is  greater  than  P^.  Path  3  represents  the  isentropic  (or 
constant-s)  compression  process  between  P^  and  P^  and  path  4 
represents  the  adiabatic  compression  between  the  same  two  pressures. 
Thus,  the  definition  of  the  adiabatic  efficiency  can  be  written  as 


n 


a 

c 


H.  ,  .  -  H, 

isentropic  28  1 


adiabatic  2  1 


Figure  6.  H-s  Diagram  for  Compression  Processes 
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Because  of  the  fact  that  isentropic  processes  cause  no  increase  in  entropy 
and  that  adiabatic  processes  do  result  in  increased  entropy,  H2  “ 

will  always  be  greater  than  or  equal  to  112^  “  for  adiabatic 

compression  processes  because  the  constant  pressure  lines  diverge  on  the 
H-s  diagram.  The  adiabatic  efficiency  of  the  processes,  therefore,  will 
always  be  less  than  one.  The  component  adiabatic  efficiencies  for  the  fan 
and  the  compressor  are  user  input  parameters  and  are  assumed  to  be 
constant  at  all  operating  points. 

Two  subroutines  are  used  to  simulate  the  effects  of  the  compression 
devices:  CMPRSS  and  CMPRSR.  CMPRSS  calculates  the  enthalpy  required  to 

perform  a  specified  compression  with  a  given  adiabatic  efficiency. 
CMPRSR,  on  the  other  hand,  calculates  the  pressure  rise  based  on  the 
enthalpy  input  to  the  compression  device  and  its  component  adiabatic 
efficiency . 

3.  The  Burner  Section:  Airflow  exiting  the  compressor  has  four 

cooling  flows  ducted  off  to  be  used  for  high  and  low  pressure  vane  and 
rotor  cooling  in  the  turbines  before  it  enters  the  burner  section  of  the 
engine  at  station  31.  When  it  reaches  station  33  it  encounters  a  splitter 
plate.  Some  of  the  flow  then  goes  through  the  stoichiometric  combustion 

chamber  from  station  33  to  36  and  the  rest  is  ducted  around  for  use  as 

dilution  air.  The  dilution  air  and  the  products  of  stoichiometric 
Combustion  are  then  mixed  to  form  the  burner  exit  flow. 
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Figure  7.  The  Burner  Section  Model 

The  Dilution  Fraction:  The  fraction  of  total  burner  air  used  for 
dilution  (PMBDL  in  the  program:  Percent  Massflow,  B^urner  £iLution)  is 
determined  in  the  following  ways  for  the  three  modes  of  engine  operation; 

a.  Design:  In  the  DESIGN  subroutine  this  fraction  is  determined 
through  an  iterative  process  so  that  the  temperature  of  the  flow  that 
results  from  the  mixing  of  dilution  air  with  the  stoichiometric  combustor 
exit  is  at  its  specified  maximum  allowable  value.  This  maximum  allowable 
value  can  be  specified  in  either  of  two  ways. 

(1)  For  a  control  scheme  based  on  gas  temperature,  the 
temperature  exiting  the  burner  at  station  38  can  be  specified  directly 
since  Che  flow  from  station  38  Co  station  40  where  it  first  enters  Che 
turbine  vanes  is  assumed  to  be  adiabatic. 
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(2)  For  a  control  scheme  based  on  turbine  metal 
temperature,  the  maximum  allowable  ^as  temperature  exiting  the  burner  is 
calculated  from  the  maximum  allowable  metal  temperature,  the  temperature 
and  amount  of  cooling  flow  being  used  for  the  first  stage  vanes  in  the 
high  pressure  turbine  and  the  cooling  effectiveness  for  the  manner  in 
which  they  are  being  cooled. 

b.  Off-Design:  For  cases  where  operation  at  an  off-design 

flight  condition  is  to  be  calculated,  the  burner  dilution  fraction  which 
causes  one  of  the  control  parameters  to  meet  its  constraint  (i.e.,  fan  or 
compressor  physical  or  corrected  speed  reaching  lOOX  of  its  design  value 
or  first  stage  turbine  vane  metal  or  turbine  inlet  gas  temperature 
reaching  its  design  limit)  while  all  the  others  are  less  than  their 
limiting  values  is  determined  through  an  iterative  routine,  OFFDSG. 

c.  Throttled  Power  Settings:  The  fraction  of  burner  air  used 
for  dilution  is  incrementally  increased  from  its  value  at  intermediate 
power  to  a  calculated  maximum  to  simulate  peiformance  at  the  requested 
number  of  throttled  power  settings.  Each  incremental  step  is  one 
twentieth  of  the  difference  between  PMBDL  at  intermediate  power  and  the 
calculated  maximum.  The  maximum  value  of  the  burner  dilution  fraction 
corresponds  to  an  "idle"  power  setting.  PMBDL  for  "idle"  is  calculated 
with  a  formula  in  subroutine  THROTL  which  was  derived  "empirically" 
through  many  runs  of  the  program  TRBOFN  for  many  different  engine  designs 
to  determine  the  maximum  burner  dilution  fraction  that  allowed  the  engines 
to  run. 

The  Stoichiometric  Combustor:  The  stoichiometric  combustion  chamber 
in  the  engine  model  is  the  section  between  stations  33  and  36.  Between 
stations  33  and  3A  heat  is  being  transferred  from  the  air  to  the  fuel 
until  their  temperatures  are  equal.  From  station  34  to  36,  the 
stoichrometric  chemical  reaction  between  the  fuel  and  the  oxygen  in  the 
air  takes  place.  The  thermodynamic  effects  of  this  reaction  are  accounted 
for  in  subroutine  STOICH.  The  stoichiometric  fuel-to-air  ratio  is 
determined  by  the  ratio  of  carbon  and  hydrogen  atoms  in  the  fuel  molecules 
to  the  number  of  oxygen  atoms  required  to  complete  the  chemical  reaction 
from  fuel  and  oxygen  in  the  air  to  water  vapor  and  carbon  dioxide.  The 
stoichoimetric  burning  of  the  fuel  results  in  no  oxygen  molecules  being 
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left  in  the  products  of  combustion  at  station  36.  In  the  burner  model  the 
specific  enthalpy  of  the  products  of  stoichiometric  combustion  is 
determined  based  on  the  fuel  lower  heating  value  and  the  stoichiometric 
fuel-to-air  ratio.  This  quantity  is  used  as  one  of  the  inputs  to  the 
calculation  in  which  the  effects  of  mixing  dilution  air  ducted  from  the 
compressor  exit  are  determined. 

Dilution  Air  Mixing:  The  mixing  of  the  dilution  air  stream  (station 
30)  with  the  hot  gas  stream  exiting  the  stoichiometric  combustion  chamber 
at  station  36  to  form  the  stream  at  station  38  may  be  performed  by  one  or 
more  members  of  a  set  of  subroutines.  When  the  desired  exit  gas 
temperature  i-s  known,  as  in  the  design  case,  then  the  subroutine 
DILUTE  iterates  to  find  the  burner  dilution  air  fraction  that  results  in 
the  exit  stream  having  the  specified  temperature.  If,  on  the  other  hand, 
the  burner  dilution  air  fraction  is  specified  (as  at  throttled  power 
settings  and  during  iterations  to  find  the  intermediate  power  setting  at 
some  off-design  flight  condition)  then  the  MIX...  subroutines  performs  the 
calculations  to  mix  the  two  streams. 

4.  The  Turbine  Section:  Both  the  high  pressure  turbine  and  the  low 
pressure  turbine  are  modeled  in  exactly  the  same  way.  Figure  8  is  the 
schematic  representation  of  the  high  pressure  turbine  model. 
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Station  40  is  the  turbine  section  entrance.  The  turbine's  vane  cooling  is 
modeled  as  a  constant  area  mixer  in  which  the  main  stream  at  station  4A  is 
mixed  with  cooling  air  which  was  ducted  around  the  burner  section  from 
just  behind  the  compressor  to  station  4B.  The  mixture  of  these  two 
streams  forms  the  stream  at  station  4C.  Station  41  is  the  choke  point  in 
the  first  stage  vanes  which  is  modeled  as  a  constriction.  Station  42  is 
the  entrance  to  the  work  extraction  section,  the  rotor,  which  ends  at 
station  40.  Station  40  is  also  the  entry  plane  of  the  mainstream  to  the 
rotor  cooling  model.  The  mixing  of  cooling  air  at  station  4E  with  the 
mainstream  at  40  is  modeled  in  exactly  the  same  manner  as  the  vane 
cooling.  The  exit  stream  formed  by  this  mixing  is  the  exit  stream  from 
the  turbine  since  the  flow  undergoes  no  further  changes  to  station  46. 

Basis  for  Turbine  Section  Model:  Engines  with  cooled  turbine  parts 
generally  make  use  of  similar  cooling  schemes.  Typical  designs  of  cooled 
turbine  guide  vanes  make  use  of  film  cooling.  The  cooled  vanes  are 
upstream  of  the  section  that  performs  the  work  extraction,  the  rotor. 
Moat  rotor  cooling  schemes  rely  on  internal  convection  to  do  most  of  the 
cooling  by  ducting  most  of  the  cooling  air  through  the  hollow  rotor  blades 
and  blowing  it  out  their  trailing  edges.  For  a  single  stage  turbine, 
then,  the  work  extraction  process  is  performed  on  gas  flow  that  includes 
the  air  that  was  used  to  cool  the  vanes  but  not  the  rotor  cooling  air. 
The  turbine  section  model  accounts  for  these  assumptions  relative  to 
realistic  cooled  turbine  designs  in  the  following  ways: 

(1)  The  constant  area  mixing  model  balances  enthalpy,  momentum 
and  mass  flow  between  the  two  input  streams  and  the  exit  stream  before  the 
work  extraction  calculation  to  simulate  the  effects  of  vane  cooling  on  the 
flow  from  which  the  rotor  must  extract  work. 

(2)  The  vane  metal  temperature  is  calculated  based  on  the 
temperature  and  amount  of  the  cooling  air,  the  temperature  of  the  hot  gas 
to  which  the  vane  metal  is  exposed,  and  the  cooling  effectiveness  of  the 
cooling  scheme.  The  cooling  effectiveness  for  each  of  the  four  cooling 
processes  in  the  turbines  are  user-selected  and  assumed  to  be  constant  for 
all  operating  points. 
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Cooling  effectiveness  is  a  measure  of  how  well  a  particular  cooling  scheme 
works  based  on  a  ratio  of  temperature  differences.  It  is  defined  as  the 
ratio  of  the  difference  between  temperatures  of  the  cooled  metal  and  the 
cooling  gas  to  the  difference  between  the  temperatures  of  the  hot  gas  and 
the  cooling  gas,  or: 


T  -  T 

hot  gas  metal 


T  -  T 

hot  gas  cool  gas 


As  an  example  of  the  correlation  between  technology  level  and  cooling 
effectiveness,  a  current  technology  engine's  high  pressure  turbine  vane 
cooling  scheme  could  result  in  a  cooling  effectiveness  of  about  73Z  using 
just  over  8%  of  the  core  airflow  for  cooling. 

(3)  After  the  work  extraction  calculations  are  performed  on  the 
stream  that  includes  the  gas  flow  entering  the  turbine  section  plus  the 
vane  cooling  flow,  the  turbine  rotor  cooling  flow  is  mixed  into  the 
stream.  This  mixing  is  modeled  by  assuming  it  takes  place  in  a  constant 
area  duct  in  which  the  mass  flow,  enthalpy  and  momentum  of  the  exit  stream 
all  balance  with  the  sums  from  the  two  entering  streams. 

Turbines  as  Work  Extraction  Devices:  Both  turbines  have  user 
selectable  component  adiabatic  efficiencies  associated  with  them.  These 
efficiencies  are  assumed  constant  regardless  of  the  operating  point  just 
as  those  associated  with  the  fan  and  compressor.  The  definition  of  the 
adiabatic  efficiency  for  an  expansion  process  is  the  ratio  of  the  enthalpy 
extracted  in  an  adiabatic  expansion  process  to  that  collected  in  an 
isentropic  expansion  with  the  same  pressure  drop.  On  Figure  9,  another 
H-s  diagram,  lines  1  and  2  are  lines  of  constant  pressure  where  is 
higher  than  P^. 
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Figure  9>  H-s  Diagram  for  Expansion  Processes 

Path  3  describes  the  isentropic  (constant-s)  expansion  from  Pj^  to  P^ 
and  path  A  describes  the  adiabatic  expansion  between  Che  same  two 
pressures.  The  definition  of  the  adiabatic  efficiency  of  an  expansion 
process  can  Chen  be  expressed  as; 

^  adiabatic  ^1  2 

*e  ~ir  I  ;  H,  -  H, 

isentropic  1  2s 

Since  the  most  efficient  expansion  process  is  the  constanC-s  path  shown  as 
path  3,  and  inefficient  processes  result  in  increased  entropy,  the  turbine 
adiabatic  efficiency  must  always  be  less  Chan  or  equal  to  one.  Two 
subroutines  perform  the  calculations  to  simulate  Che  effects  of  Che 
turbines:  TURBIN  and  TURBND.  TURBND  is  used  to  calculate  Che  turbine 
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expansion  ratio  when  the  work  output  required  (delta  enthalpy,  H)  is  known 
and  the  adiabatic  efficiency  is  specified.  This  is  used  in  the  design 
case  when  the  compressor  and  fan  pressure  ratios  are  specified.  TURBIN 
calculates  the  work  output  based  on  a  given  value  for  the  turbine 
expansion  ratio  and  a  specified  component  adiabatic  efficiency.  This  is 
used  for  all  off-design  operating  conditions. 

Turbine  Vane  and  Rotor  Cooling;  All  flow  cooling  processes  in  the 
turbines  are  modeled  in  exactly  the  same  way  by  the  MIX...  subroutines. 

Members  of  that  series  of  subroutines  calculate  total  and  static 

thermodynamic  properties  of  the  mixture  of  two  streams  based  on: 

1.  Known  totals  and  static  conditions  and  the  constituent 

make-up  of  the  main  stream. 

2.  Known  total  conditions  in  the  cooling  stream. 

3.  Known  flow  fraction  used  for  cooling  (assumed  constant) 

4.  Requiring  the  static  pressures  of  the  two  streams  to  be  mixed 
to  be  the  same  at  the  mixing  plane. 

5.  Conservation  of  Massflow. 

6.  Conservation  of  Energy  (to  get  the  total  conditions). 

7.  Force  -  Momentum  Balance  (to  get  the  static  conditions). 

Force  ^  0  since  walls  are  parallel  due  to  the  constant  area  assumption. 

5.  The  Augmentor:  Only  mixed  flow  engines  can  be  augmented  in  this 
program.  The  augmentor  model  is  very  similar  to  the  burner  model  in  that 
the  afterburner  of  the  engine  is  treated  as  two  ducts:  one  for 
stoichiomertic  burning  and  one  for  dilution.  The  major  difference  is  that 
the  input  stream  to  the  afterburner  is  not  pure  air  but  rather  a  stream 
that  contains  some  products  of  combustion  already.  Figure  10  is  a 
schematic  representation  of  the  afterburner  model.  Note  the  similarity 
between  this  and  Figure  7,  the  burner  section  model. 
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15  6D  6E  8B 


Figure  10.  The  Afterburner  Model 


Stations  63  through  66  are  analogous  to  stations  33  through  36  in  the 

burner  model.  The  fuel  is  heated  by  the  air  between  stations  63  and  64 

and  it  is  burned  stoichiometrically  with  the  oxygen  in  the  entry  stream 
between  stations  64  and  66.  Subroutine  AFBKMR  calls  subroutine  STOICH  to 
perform  these  calculations.  The  MIX...  subroutines  are  used  to  account 
for  the  mixing  of  the  stoichiometric  combustor  exit  and  the  dilution 
streams.  Likewise,  these  MIX...  routines  are  used  to  account  for  the 
return  of  the  nozzle  liner  and  flap  cooling  streams  at  stations  6E  and  SB. 

6.  The  Nozzle:  The  nozzle  is  modeled  as  a  reduced  area  section  for  the 

flow  to  pass  through.  Its  throat  area  is  calculated  in  subroutine  DESIGN 

as  that  which  chokes  the  flow  at  the  design  flight  condition.  At 
off-design  flight  conditions  it  is  not  always  choked  but  it  does,  however, 
always  adjust  the  flow  for  perfect  expansion  into  the  atmosphere  at  all 
flight  condition,  including  design.  In  other  words,  the  throat  area,  A80, 
is  always  held  constant  for  throttled  or  off-design  operation  but  the  exit 
area,  A90,  varies  appropriately  so  the  flow  is  always  perfectly  expanded. 
This  also  holds  true  for  the  bypass  nozzle  for  separate  flow  engines. 
When  the  afterburner  is  turned  on,  both  the  nozzle  throat  area  and  the 
exit  area  are  changed  to  allow  the  additional  massflow  to  pass  and  to  be 
perfectly  expanded.  The  MIX...  subroutines  are  used  to  model  the  return 
of  the  nozzle  liner  and  flap  cooling  flows  in  mixed  flow,  augmented 
engines.  When  these  nozzle  cooling  ducts  exist  they  are  assumed  to  always 
remain  open  for  passing  flow  whether  or  not  the  augmentor  is  turned  on. 
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DEFINITIONS  OF  ENGINE  EFFICIENCIES 


DEFINITIONS  OF  ENGINE  EFFICIENCIES:  Some  of  Che  performance  outputs  from 
TRBOFN  are  engine  efficiencies  and  thrust  specific  fuel  consumption. 
Their  definitions  and  methods  of  calculation  follow. 

1.  TSFC  Definition:  Thrust  Specific  Fuel  Consumption  (TSFC)  is 
defined  as  the  rate  at  which  fuel  is  being  burned  divided  by  the  amount  of 
thrust  being  produced  as  a  result.  That  is: 

TSFC  =  m,/F 
t  n 

The  units  on  TSFC  are  pounds  of  fuel  per  hour  per  pound  of  thrust. 

2.  Efficiency  Definitions:  Defining  the  engine  efficiencies 
described  below  will  help  show  the  dependence  of  TSFC  on  how  efficiently 
the  engine  uses  the  energy  released  by  the  burning  of  fuel. 

Thermal  Efficiency;  Thermal  efficiency  is  a  measure  of  how 
well  the  core  of  the  engine  makes  use  of  the  energy  made  available  by 
burning  fuel.  Only  some  of  the  air  acted  upon  by  the  fan  goes  through  the 
core  of  the  engine.  The  rest  goes  through  the  bypass  duct  (Figure  11). 
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Figure  11.  Part  of  the  Low  Pressure  Turbine  Drives  the 
Part  of  the  Fan  that  Acts  on  the  Bypass  Air 


The  work  required  to  drive  only  the  portion  of  the  fan  that  acts  on  the 
air  that  goes  through  the  core  of  Che  engine  is  a  fraction  of  the  total 
work  extracted  by  the  low  pressure  turbine.  We  can  think  of  this  as 
being  done  by  a  turbine  Chat  is  only  a  fraction  as  large  as  the  original 
turbine  which  was  required  Co  extract  a  sufficient  amount  of  work  to 
drive  Che  entire  fan.  Consider  Che  engine  with  these  two  "partial" 
components  in  it  (Figure  12). 


Figure  12.  Partial  Components  for  Thermal  Efficiency  Calculation 

Thermal  efficiency  is  defined  as  Che  amount:  of  energy  available  Co  do 

useful  work  that  remains  in  the  core  airstream  after  a  sufficient  amount 
of  work  has  been  drawn  from  the  stream  by  the  portion  of  the  low  pressure 

turbine  needed  to  drive  the  part  of  the  fan  which  acts  only  on  the  air 

Chat  goes  through  Che  core  of  Che  engine.  The  way  Co  calculate  the  rate 
at  which  energy  is  made  available  Co  do  useful  work  is  Co  calculate  Che 
rate  at  which  energy  can  be  extracted  from  the  stream  with  an  ideal 
turbine,  i.e.  a  turbine  with  a  component  efficiency  of  unity  while 
leaving  a  sufficient  amount  of  energy  to  produce  zero  net  thrust  on  the 
system  after  perfect  expansion  through  a  nozzle.  This  device  is 

represented  in  Figure  12  as  Che  large  cross-hatched  turbine  labeled  ITl 
(Ideal  Turbine  1).  The  zero  thrust  requirement  in  Che  definition  makes 
thermal  efficiency  independent  of  Che  flight  velocity.  The  rate  at  which 
energy  is  released  by  Che  burning  fuel  is  calculated  as  Che  product  of 
the  mass  flow  of  Che  fuel  and  the  fuel  lower  heating  value.  The  fuel 
lower  heating  value  is  the  amount  of  energy  released  per  mass  of  fuel 
burned.  Its  units  are  Btu  per  pound  mass. 
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Transfer  Efficiency:  Transfer  efficiency  is  a  measure  of 
how  efficiently  energy  is  transferred  from  the  low  pressure  turbine  to 
the  bypass  stream.  It  is  defined  as  the  total  energy  in  both  the  core 
and  bypass  streams  available  to  do  useful  work  after  the  low  pressure 
turbine  has  drawn  off  enough  energy  to  drive  the  entire  fan  divided  by 
the  work  drawn  off  by  the  ideal  turbine  described  in  the  thermal 
efficiency  definition.  Referring  to  Figure  13,  the  numerator  is 
calculated  as  the  sum  of  the  work  extracted  by  two  ideal  turbines;  one  in 
the  core  stream,  IT2,  and  the  other  in  the  bypass  stream,  IT3.  The 
energy  extracted  by  each  of  these  turbines  is  determined  such  that  the 
exit  streams  each  produce  zero  net  thrust  after  perfect  expansion  through 
a  nozzle.  The  denominator  is  the  work  extracted  by  ITl  in  Figure  12. 
That  is:  ^^2  ^  ^3 


Figure  13.  Ideal  Work  Extraction  for  Transfer  Efficiency  Calculation 
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Propulsive  Efficiency:  Propulsive  efficiency  is  a  measure 
of  how  well  energy  available  is  converted  to  thrust  power.  It  is 
calculated  as  the  fraction  of  the  total  energy  available  to  do  useful 
work  in  both  exhaust  streams  that  shows  up  as  actual  thrust  work.  The 
total  energy  in  the  exhaust  streams  is  the  sum  of  the  work  extracted  by 
the  two  ideal  turbines  used  in  the  transfer  efficiency  definition,  IT2 
and  IT3  in  Figure  13.  Thrust  power  is  simply  the  thrust  (force)  times 
the  flight  velocity.  That  is: 


F  u 
n  o 


pr 


W2+W3 
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THE  DESIGN  PROCESS 
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THE  DESIGN  PROCESS:  You  should  find  it  easy  to  follow  the  design  process 
for  both  mixed  and  separate  flow  engines  in  detail  by  reading  the 
comments  and  FORTRAN  code  that  make  up  subroutine  DESIGN.  The  preceding 
explanations  of  the  basis  and  workings  of  the  component  models  should 
make  each  of  them  very  easy  to  understand  when  reading  the  code  and 
comments  that  make  them  up. 


75 


OFF-DESIGN  OPERATION:  Off-design  operation  of  a  turbofan  engine  includes 
throttled  power  settings  and  intermediate  power  settings  at  other-than- 
design  flight  conditions.  Calculating  off-design  performance  for  a 
turbofan  engine  is  considerably  more  complicated  than  calculating  design 
performance  because  balanced  off-design  engine  performance  can't  be 
calculated  analytically  with  the  assumptions  we  used  to  generate  this 
model.  We  chose  to  use  a  "nested  DO  loop"  approach  in  iterating  to  find 
balanced  off-design  operating  points.  Several  other  approaches  have  been 
used  in  other  engine  performance  simulation  programs,  most  of  which 
require  less  programming  steps  and  run  faster,  but  we  have  found  them 
terribly  difficult  to  understand  with  the  poor  documentation  available 
with  them.  While  the  method  of  programming  we  chose  requires  more  coding 
and  results  in  longer  run  times,  it  also  makes  it  possible  to  write  an 
understandable  program  that  you  don't  have  to  be  the  author  of  to 
understand  and  use.  We  authors  were  willing  to  sacrifice  running  speed 
to  write  an  easy-to-use  program  that  users  could  understand.  We  followed 
a  comprehensive  set  of  documentation  standards  in  writing  every  part  of 
the  system  of  which  TRBOFN  is  a  part.  Because  of  the  extensive 
documentation  in  the  programming  itself,  we  believe  that  all  the 
subroutines  that  make  up  this  system  should  be  easy  to  understand  just  by 
reading  the  "Purpose  statement"  and  the  comment  lines  throughout  the 
coding. 

1.  Throttling:  Balanced  throttled  power  settings  are  found  by  the 
subroutine  THROTL.  THROTL  incrementally  increases  the  burner  dilution 
fraction  to  simulate  throttled  power  settings.  It  then  calls  TRBOFN's 

workhorse  subroutine,  BALENG,  which  gets  its  name  from  what  it  does: 

BALance  the  ENGine.  Balancing  the  engine  refers  to  finding  the  set  of 
engine  cycle  parameters  that  satisfy  conservation  of  energy,  momentum  and 
massflow  throughout  the  engine.  BALENG  sets  initial  guesses  on  several 
parameters  it  needs  to  get  started  and  calls  three  other  major 

controlling  subroutines  which  are  described  below.  You  can  see  which 

parts  of  the  engine  each  of  these  major  control  routines  operates  on  by 
looking  at  Figure  14. 
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Figure  14b.  Parts  of  Separate  Flow  Engines  Controlled  by  the 
Major  Subroutines  Called  by  BALENG 


THROTl  performs  calculations  using  appropriate  component  models 
to  ensure  balanced  operation  of  the  engine  between  the  compressor  exit, 
station  31,  and  the  choke  point  in  the  high  pressure  turbine’s  first 
stage  vanes,  station  41. 

THR0T2  iterates  to  find  a  balance  between  stations  41  and  47, 
the  choke  points  in  the  high  and  low  pressure  turbines'  first  stage 
vanes.  One  resultant  parameter  from  THR0T2  is  the  work  extracted  from 
the  flow  by  the  high  pressure  turbine.  BALENG  uses  this  to  iterate  to 
balance  the  core  of  the  engine.  That  is,  BALENG  changes  the  driving 
parameter  in  its  core  balance  loop  until  the  work  extracted  by  the  high 
pressure  turbine  balances  with  the  work  required  to  drive  the  compressor 
at  the  current  compressor  pressure  ratio. 

THR0T3  is  called  each  time  the  core  balance  loop  is  satisfied. 
Its  job  is  to  calculate  the  work  extracted  by  the  low  pressure  turbine 
and  to  do  either  of  two  other  things,  based  on  whether  the  engine  design 
is  mixed  flow  or  separate: 

For  mixed  flow  designs:  THR0T3  iterates  to  find  the 
appropriate  low  pressure  turbine  expansion  ratio  which  results  in  the 
required  static  pressure  balance  for  mixing  the  core  exit  stream  at 
station  50  and  the  bypass  stream  at  station  15. 

For  separate  flow  engines:  THR0T3  iterates  to  find  the  low 
pressure  turbine  expansion  ratio  which  results  in  a  choked  core  nozzle  if 
the  back  pressure  exerted  by  the  atmosphere  is  sufficiently  low  to  allow 
choking.  If  the  pressure  in  the  core  duct  relative  to  that  in  the 
atmosphere  is  insufficient  to  choke  the  nozzle,  then  THR0T3  iterates  to 
find  the  appropriate  turbine  expansion  ratio  that  results  in  a  perfectly 

expanded  nozzle  with  exit  area  equal  to  the  design  choke  area,  A„„. 
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In  either  case,  THR0T3  returns  the  work  extracted  by  the  low  pressure 
turbine  to  BALENG  which  uses  this  quantity  to  compare  to  the  work 
required  to  drive  the  fan  and  adjusts  the  fan  pressure  ratio  to  match  the 
work  extracted  by  the  turbine  to  get  the  fan  pressure  ratio  for  the  next 
iteration.  BALENG  continues  to  iterate  in  this  manner,  each  time  through 
iterating  to  achieve  a  balanced  core  and  adjusting  fan  pressure  ratio 
based  on  the  information  THR0T3  returns  until  the  engine  is  balanced  all 
the  way  through  (i.e.  the  work  out  from  the  high  pressure  turbine  matches 
the  work  required  to  run  the  compressor,  an  analogous  situation  is  true 
for  the  fan  spool,  the  first  stage  vanes  of  both  turbines  are  choked  and 
have  the  appropriate  massflows  going  througii  them  and  the  nozzle  is 
choked  and/or  perfectly  expanded). 

2.  Off-Design  Flight  Conditions:  To  find  the  intermediate  power 
setting  at  an  off-design  flight  condition,  XRBOFN  calls  subroutine  OFFDSG 
which  calls  BALENG  to  performs  Che  functions  described  above  under 
Throttling.  Initial  guesses  for  the  iterative  process  that  leads  to 
finding  Che  off-design  flight  condition’s  intermediate  power  setting  are 
generated  in  subroutine  INTGES  (INiTial  GuESs).  OFFDSG  (OFF  DeSiGn)  uses 
subroutine  BEST05  (BEST  Of  5)  to  get  consecutive  guesses  on  the  burner 
dilution  fraction  (PMBDL)  in  its  iteration  to  find  the  intermediate  power 
setting.  Once  the  intermediate  power  setting's  PMBDL  is  known,  finding 
throttled  power  settings  at  off-design  flight  conditions  is  done  in 
exactly  Che  same  way  as  at  Che  design  flight  condition. 

a.  The  Intermediate  Power  Setting:  The  intermediate  power 
setting  at  an  off-design  flight  condition  is  defined  as  the  state  of 
balanced  engine  operation  where  one  of  Che  possible  control  parameters  is 
at  its  control  limit  while  all  the  others  are  below  theirs. 
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b.  The  Possible  Control  Limits:  There  are  five  possible 
control  limits  which  must  be  checked  for  each  balanced  state  of  engine 
operation  (i.e.,  each  balanced  engine  operating  condition  found  by  BALENG 
for  each  guess  on  PMBDL  generated  by  BEST03  from  OFFDSG) .  The  five 
possible  control  limits  are: 

1)  The  physical  rotational  speed  of  the  low  pressure  (fan) 
spool  (PPNl)  is  calculated  as  a  fraction  of  the  design  speed  by 
subroutine  OPPCN  (Operating  Point  Corrected  N,  where  N  is  the  spool  speed 
in  RPH).  Its  limiting  value  is  lOOZ  of  the  design  speed. 

2)  The  physical  rotational  speed  of  the  high  pressure 
(compressor)  spool  (PPN2)  is  calculated  and  treated  exactly  like  the  fan 
physical  spool  speed. 

3)  The  corrected  rotational  speed  of  the  low  pressure 
spool  (PCNl)  is  also  calculated  as  a  fraction  of  the  design  value  by 
subroutine  OPPCN.  It  too  is  limited  to  1002  of  the  design  value. 

4)  The  corrected  rotational  speed  of  the  high  pressure 
spool  (PCN2)  is  calculated  and  treated  exactly  like  the  corrected  speed 
of  the  fan  spool. 

5)  The  fifth  possible  control  parameter  is  selected  by  the 
user  during  the  design  process  as  either: 

a)  T40G,  the  gas  temperature  entering  the  high 
pressure  turbine's  first  stage  vanes  (Turbine  Inlet  Temperature), 

b)  or  T40MVA,  the  high  pressure  turbine's  first  stage 
Vane  Metal  temperature  Allowable. 

The  user  also  sets  the  design  (limiting)  value  of  this  parameter  during 
the  design  process. 
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c)  BESX05:  Each  time  OFFDSG  calls  BEST05«  it  projects  a  new 
guess  on  PMBDL  for  each  of  the  five  control  parameters  by  assuming  the 
variation  of  each  control  parameter  is  linear  in  PMBDL.  These  five 

linear  projections  result  in  five  different  values  of  PMBDL,  each  of 

which  would  put  one  of  the  control  parameters  at  its  limiting  value  if 
the  parameter's  value  really  varied  linearly  with  PMBDL.  The 

intermediate  power  setting  is  defined  as  the  point  of  balanced  engine 

operation  where  only  one  of  the  control  parameters  is  ^  its  limit  while 
all  the  others  are  BELOW  theirs.  BEST05,  therefore,  returns  only  the 
highest  value  of  PMBDL  of  the  five  it  calculates  to  be  used  as  the  next 
guess  to  be  used  in  BALENG.  It  chooses  the  highest  value  of  PMBDL 

because  PMBDL  is  the  percent  of  the  total  burner  air  used  for  dilution 
and  a  higher  value  of  PMBDL  means  less  fuel  is  being  burned.  When  more 
fuel  is  burned,  the  engine  runs  hotter  and  faster.  Setting  PMBDL  at  one 

of  the  lower  values  calculated  by  BEST05  would  result  in  at  least  one 

control  parameter  having  a  value  larger  than  its  limiting  value  and 

therefore  would  not  correspond  to  the  intermediate  power  setting. 
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SUBROUT ISE  GROUPS 
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OTPi 


ALPHABETICAL  LIST  OF  SUBROUTINES  BY  GROUPS : 


Gas  Model 


Comp  Model  Control 


Utility  Performance 


AFBRNR  X 

ATMOS  X 

BALENG 

BEST05 

CMP...  X 

DESIGN 

DILUTE  X 

ETA 

GASPAC 

LOADGl  X 

MFRACS  X 

THERMH  X 

THERM!  X 

THERMP  X 


IDLXIT 

INPUT 

INTGES 

LOADAl 

MIX...  X 

NEWGES 

N02C0N  X 

OFFDSG 

OPPCN 

OUT... 

PDROP  X 

RAM 

RAMREC  X 

SETEQ 

STAT... 

STOICH  X 

THROT. . . 

THRUST 

TRBOFN* 

TURB...  X 
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X 

X 

X 

X 

X 

X 

X 

X 


*  TRBOFN  is 


the  name  of  the  mainline  program 
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Chapter  4. 

ENCOM:  The  Engine  Sizing  and  Comparison  Program 
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1.  Introduction:  ENCOM  is  an  interactive  program  using  a  question-answer 
format  to  get  input  information  from  the  user.  It  can  read  engine 
performance  data  from  up  to  five  files  generated  by  the  program  TRflOFN  on 
each  run.  That  means  it  can  be  used  to  compare  up  to  five  engines  at  a 
time  or  to  generate  a  file  to  plot  a  trend  with  five  points. 

ENCOM  has  three  basic  functions: 

1.  SIZING:  to  determine  the  minimum  size  of  each  engine  relative  to 
the  others  for  each  to  be  able  to  meet  or  exceed  each  of  up  to 
five  user  input  maneuver  requirements, 

2.  COLLECTING  DATA  FOR  PLOTS:  to  collect  the  requested  data  from 
each  of  the  five  engine  performance  data  files  so  the  user  can 
make  plots  of  any  stored  parameter  versus  any  other,  and 

3.  MISSION  ANALYSIS:  to  calculate  performance  throughout  and  to 
ultimately  determine  the  minimum  fuel-to-take  off  gross  weight 
fraction  necessary  for  each  engine  to  complete  an  eleven  leg 
fighter/ground  attack  mission  in  the  same  fighter  aircraft. 

The  purpose  of  this  chapter  is  to  tell  you  about  what  ENCOM  can  do.  To 
find  out  exactly  how  ENCOM  does  all  these  things  you'll  have  to  read  the 
documentation  included  in  the  actual  FORTRAN  source  coding  of  the 
program.  Like  all  the  other  programs  in  this  system,  ENCOM  is  made  up  of 
a  mainline  control  code  which  calls  several  subroutines  to  perform 
specific  functions.  The  mainline  and  every  subroutine  follow  the  same 
documentation  standards:  each  includes  a  statement  of  purpose  followed  by 
an  explanation  of  each  variable  passed  into  and  out  of  the  subroutine  and 
at  least  one  comment  line  for  every  line  of  FORTRAN  coding  throughout  the 
program. 
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SIZING: 

The  sizing  process  in  ENCOM  requires  several  user  inputs; 

o  up  to  five  maneuver  requirements  (i.e.,  flight  Mach  number, 
altitude,  sustained  acceleration  requirement,  take-off  wing  loading  and 
maneuver  wing  loading) 

o  which  engine  is  to  be  used  as  the  standard  for  comparison, 

o  the  names  of  the  data  files  created  by  TRBOFN  to  be  read  by  ENCOM 
to  get  the  engine  performance  information  needed  to  size  the  engines 
relative  to  each  other.  NOTE:  All  engine  performance  data  files  to  be 
used  by  ENCOM  must  have  sea  level  static  performance  information  . 

Once  one  engine  is  selected  as  the  standard  for  comparison,  the  other 
engines'  sizes  are  scaled  up  or  down  relative  to  that  engine's  so  that 
their  thrust  output  at  each  maneuver  requirement  is  sufficient  to  propel 
the  aircraft  through  Che  specified  maneuver.  This  results  in  a  sizing 
factor  beign  calculated  for  each  engine  relative  to  the  standard  for  each 
maneuver  requirement.  It  also  enables  ENCOM  to  calculate  the  sea  level 
static  thrus t-to-we ighc  ratio  for  the  airplane  with  each  engine  scaled  by 
each  sizing  factor.  The  critical  maneuver  requirement  is  the  one  that 
forces  the  airplane  to  have  Che  tiighest  sea  level  static  thrus t-to-we ight 
ratio.  That  maneuver  requires  the  engine  to  be  scaled  up  the  most.  If 
the  engine  is  scaled  to  meet  the  critical  maneuver  requirement,  then  it 
will  be  more  than  big  enough  to  complete  the  other  maneuver  requirements. 
Once  the  critical  maneuver  requirement  is  determined  for  each  engine,  the 
corresponding  sizing  factor  is  used  to  calculate  that  engine's  physical 
size,  total  airflow  size  and  core  airflow  size  relative  to  the  standard's. 
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COLLECTING  DATA  FOR  PLOTS; 


Subroutine  PLOTCON  is  called  by  ENCOM  to  control  the  collection  of  data 
for  plots.  PLOTCON  uses  other  subroutines  to  perform  parts  of  the 
collection  functions.  For  instance,  WPATIPLOT  is  the  subroutine  that 
displays  the  choices  for  parameters  in  the  engine  performance  data  files 
that  PLOTCON  can  get  and  output  in  the  proper  form  for  use  by  GRAPH. 
WHAT2PL0T  also  prompts  the  user  to  input  his  choicec  for  the  parameter  he 
wants  on  each  exis  of  the  graph.  The  parameters  are  displayed  as  a  table 
of  rows  and  columns  and  the  user  specifies  which  parameter  he  wants  by 
entering  its  row  and  column  location  on  the  table.  PLOTCON  then  uses 
other  subroutines  to  get  the  user  to  tell  it  whether  he  wants  throttling 
trends  plotted  for  each  engine  or  ]ust  intermediate  power  points, 
augmented  performance,  or  dry  power  only,  one  flight  condition  or  a 
series,  etc.  Once  PLOTCON  know^  enough  about  what  the  user  wants  it  to 
get  from  the  engine  performance  data  files,  it  calls  COLLECT... 
subroutines  to  do  the  data  collection.  The  COLLECT...  subroutines  return 
the  data  requested  by  the  user  to  PLOTCON  in  arrays  which  are  then  written 
to  a  user-named  output  file  by  subroutine  WRITOUT. 
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UNCLASSIFIED 


Mission  AMALYSIS; 


ENCOM's  last  function  is  calculating  and  coaparing  the  perforaance  of  up 
to  five  engines  while  propelling  an  airplane  throu^  an  eleven  leg 
fighter/ground  attack  aission.  A  acheaatic  of  the  aiaaion  is  shown  in 
Figure  15  with  a  description  of  each  leg  of  the  aission.  A  saaple  aission 
output  file  also  follows. 


Alt  (ft) 


Figure  1.  Schematic  of  Eleven  Leg  Fighter/Ground  Attack  Mission 


LEG  Altitude  (ft)  Flight  Mach  Distance  (n  ai)  Time  (min) 


1  Warm  Up 

SL 

static 

- 

15 

2  Takeoff 

SL 

calc  M 

to 

calc 

calc 

3  Climb  &  Accel 

SL  to  30K 

0  to  .88 

calc 

calc 

4  Outbound  Cruise 

30K 

.88 

275 

32 

5  Loiter 

25K 

.54 

sust  1.5  g  turns 

10 

6  Dash  &  Penetrat 

lOK 

1.10 

150 

13 

7  Combat 

2R 

.93 

2  360  deg  turns 

calc  for  3g 

8  Return  Dash 

2K 

1.12 

100 

8 

9  Climb  &  Decel 

2K  to  30R 

1.12  to  .88 

calc 

calc 

10  Inbound  Cruise 

30K 

.88 

325 

37 

11  Landing  Reserve 

IK 

.19 

- 

20 
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De«i£n  Point  Perforaance  Compg‘iaon 

In  order  to  verify  the  accuracy  of  the  performance  calculated  by  TRBOm, 
we  decided  to  compare  ita  outputs  to  an  accepted  standard  turbofan  engine 
performance  simulation  program  currently  in  use  at  the  Aero-Propulsion  Lab 
at  Hright-Patterson  Air  Force  Base,  Ohio:  SMOTE. 

This  section  of  the  Appendix  shows  how  the  performance  calculated  by 
TRBQFR  compares  to  that  calculated  by  SMOTE  at  the  design  point  for  over 
200  engine  designs.  We  chose  to  show  the  comparison  between  net  thrust 
and  thrust  specific  fuel  consumption  (TSFC)  at  the  design  point  for  each 
of  these  mixed  flow  turbofan  engines. 

Each  graph  in  this  section  shows  either  thrust  or  TSFC  information  for 
each  of  six  turbofan  engine  designs,  all  with  the  same  fan  pressure  ratio 
(FPR)  and  turbine  inlet  temperature  (T40),  and  with  compressor  pressure 
ratio  (CPE)  between  4.5  and  9.0.  The  designs  compared  include  T40s  of 
2400  R,  2600  R,  and  2800  R.  Four  different  design  flight  conditions  are 
included  in  this  section:  sea  level  static,  Mach  0.7  at  15,000  ft 
altitude,  Mach  0.9  at  15,000  ft,  and  Mach  0.9  at  the  tropopause,  36,089 
ft.  All  the  graphs  show  better  than  4X  agreement  for  every  engine. 
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Off-Design  Performance  Comparison 

This  section  of  Che  Appendix  %fes  intended  to  show  how  the  throttled  and 
off-design  flight  condition  perforsunce  calculated  by  TRBOFM  compares  to 
Chat  calculated  by  SHOTS.  He  wanted  to  show  the  comparison  between  net 
thrust  and  thrust  specific  fuel  consumption  (TSFC)  at  several  off-design 
power  points  for  each  of  several  mixed  flow  Curbofan  engines.  He 
encountered  an  insurmountable  problem  in  undertaking  this  task.  He  could 
not  get  SMOTE  to  converge  at  any  of  the  off-design  power  points  we  Cried 
to  run.  SMOTE  would  only  converge  on  design.  The  engines  we  tried  Co  run 
were  simplest  case  mixed  flow  curbofans  in  Chat  they  included  no 
proviaions  for  turbine  cooling  or  overboard  bleeds  and  after  two  months  of 
trying,  we  could  not  get  SMOTE  to  converge  at  even  one  of  Che  off-design 
points  we  Cried  to  run. 

Because  of  the  complications  we  encountered,  each  graph  in  this  section 
shows  thrust  versus  TSFC  throttling  information  for  each  of  three  flight 
conditions  calculated  by  TRBOFM  (the  design  flight  condition  and  two 
off-design  flight  conditions)  but  only  Che  design  point  as  calculated  by 
SMOTE.  The  designs  run  include  Che  same  three  turbine  inlet  temperatures 
(TAO)  as  run  in  the  design  point  comparison:  2400  R,  2600  R,  and  2800  R. 
Two  compressor  pressure  ratios  were  used  in  the  designs:  6.0  and  9.0. 
Four  different  design  flight  conditions  are  included  in  this  section:  sea 
level  static,  Mach  0.7  at  15,000  ft  altitude,  Mach  0.9  at  15,000  ft,  and 
Mach  0.9  at  the  tropopause,  36,089  ft. 
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